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ABSTRACT 

We present VLT/X-shooter spectra of 13 galaxies in the redshift range 1 < z < 6, which 
are strongly lensed by massive galaxy clusters. Spectroscopic redshifts are measured for nine 
galaxies, while three sources have redshifts determined from continuum breaks in their spec- 
tra. The stellar masses of the galaxies span four orders of magnitude between 10 7 and 10 11 
M Q and have luminosities at 1500 A rest-frame between 0.004 and 9 L* after correcting for 
the magnification. This allows us to probe a variety of galaxy types from young, low-mass 
starburst galaxies to massive evolved galaxies. The lensed galaxies with stellar masses less 
than 10 10 M Q have a large scatter compared to the fundamental relation between stellar mass, 
star formation rates and oxygen abundances. We provide a modified fit to the fundamental 
relation for low-mass, low-metallicity galaxies with a weaker dependence of the metallic - 
ity on either the star formation rate or stellar mass compared to low-redshift, high-mass and 
high-metallicity SDSS galaxies. 

Key words: galaxies: high-redshift - galaxies: distances and redshifts - galaxies: evolution 
- galaxies: abundances - Physical data and processing: gravitational lensing 



1 INTRODUCTION 

Strong gravitational lensing is an important tool for characterising 
the dark matter distribution and formation of large mass concentra- 
tions in a ACDM framework. Since massive galaxy clusters give 
rise to magnification factors of 10-50 of background galaxies, as 
well as stretching the source images into extended arcs, lensing can 
be used to characterise the physical properties of intrinsically sub- 
luminous galaxies and investigate spatially resolved properties of 
high-redshift galaxies dSwinbank et alj|2007t) . Grav itational lens- 
ing is also valuable for locating z > 6 gala xies dBradlev et al.l 
2008; iRichard et al J|20 lTbl ; Izitrin et alj|2012al) and investigate in- 
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dividual galaxies which may be representative of those responsible 
for reionisation. 

The Lyman break technique by construction selects galaxies 
with simi lar characteristics, such that their stacked spectra can be 
analysed (Shapley et al. 2003). Individual galaxies are too faint for 
spectroscopic studies with reasonable use of telescope ti me, apart 
from unusually bright Lyman break galaxies (LBGs) ferb et all 
l201fj|) . Other selection techniques have been used to study the 
fainter end of the high redshift galaxy lum i nosity function, incl ud- 
ing Lya emitters (e.g. iFvnbo et al.l 120 03 : Nilsson et al. 201lh or 
host galaxies of gamma-ray bursts (e.g. iHiorth et alj|2012l) . How- 
ever, they fail to reach the same level of detail because of the large 
amount of telescope time needed. Our current understanding of the 
physical properties of high-redshift galaxies is therefore limited to 
a small range in absolute magnitude. 

Gravitational lensing allows us to study intrinsically fainter 
galaxies with unusual spectral characteristics at any given redshift, 
although most lensed galaxies studied in detail to date still belong 
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to the high-mass, high-luminosity end. A few cases have provided 
a great deal of insight into the physical properties of individual 
z = 2-3 galaxie s, with the prime example being the lensed galaxy 
MS 1 512-cB58 dPettini et al.l200d : lTeplitz et alj2000l : |Pettini et all 
|2002|) . The past few years have seen a steady increase in the 
number of spectroscopic observations of lensed galaxies at 2 ~ 
2 either detected se r endipitously or in l a rge dedicated surveys 
dFosburv et alj|2003l; ICabanac et alj|2005l: ISwinbank et alj|2007l; 
Allam et al.l 120071: IStark et alj 12008: llin et alj 120091; iDiehl et"atl 
2009IJ 
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20 id iRigbv et alj|201ll ; IRichard et alj l201lj ; 



Wuvts et alj|2012bl)~ Whereas some z = 2 lensed galaxies are eas- 



ily investigated with low-resolution spectra from 8-m class tele- 
scopes, few are sufficiently bright for medium- or high-resolution 
spectroscopy, which is necessary for examining the interstellar 
medium of the galaxy. Rest-frame UV absorption line spectra al- 
low a detailed analysis of the interstellar me dium of galaxies with 
intrinsically very high star formation rates ( Cabanac et al. 2008 ; 
lOuider et alj2009ll2010l ; lDessauges-Zavadskv et alj|2010h . 

As the redshift around z ~ 2 represent the era of the peak 
of the star formation, the lensed galaxies are relevant to compare 
to scaling relations fo und in other galaxy samples, such as mass- 
metallicity relations dTremonti et alj|2004f) . Recently, a fundamen- 
tal relation between the mass, metallicity and SFR has been found. 
The relation has been calibrated for metal-rich (oxygen abundances 
12+log(0/H)>8.2), and relativ ely massive galaxies (ste llar masses 
above 10 9 ' 2 M Q ) at z < 0.1 dLara-Lopez et alj|2010h - At higher 
redshifts up to 2 ~ 2.5 and fo r galaxies with lower m etallicities the 
relation appears to still hold dMannucci et alj|2010h . Strong grav- 
itational lensing of intrinsically fainter and less massive galaxies 
allows us to explo re the fundam e ntal rel ation at even lowe r masses 
and metallicities. IRichard et ail d2011al) and IWuvts et al] d2012bl) 
analysed the rest-frame optical emission lines of lensed galaxies at 
2 = 1 — 3 with intrinsic low luminosities (0.1 L* ) and showed that 
the galaxies are more metal-rich than predicted from extrapolations 
of mass-metallicity relations. 

In this paper, we present shallow spectroscopic observations 
from VLT/X-shooter of 13 lensed galaxies selected towards mas- 
sive galaxy clusters at 2 ~ 0.4. Section [2] presents the data and 
Section [3] notes on each target and their redshift determination. In 
Section [4] we explore the galaxy properties from their continuum 
while in Section [5] we derive physical properties based on their 
emission lines. In Section [6] we present an extension of the cal- 
ibration at the low-mass, low-metallicity end of the fundamental 
relation for star forming galaxies. Section|7]presents a summary. 



2 OBSERVATIONAL DATA 

One of the science goals for the 2 nd generation VLT instrument X- 
shooter is to determine redshifts for faint single targets and analyse 
emission line galaxies at 2 > 1. X-shooter is a multi-wavelength 
cross-dispersed echelle spectrograph which splits up the light by 
dichroics into three arms: UV-blue (UVB), which covers 3000- 
5600 A, visible-R (VIS), which covers 5600 -10000 A, and near-IR 
(NIR ), which co vers from 1 to 2.5 micron (D'O dorico et al 1 l2006l : 
Vern et et al ■1201 lh . The slit length of X-shooter is 11" and the cho- 
sen slit widths determine the spectral resolution. For the observa- 
tions the slit widths were set to match the optical seeing measured 
at the telescope as listed in Table Q] Slit widths of 1.0/0.9/0.9 arc- 
sec in the UVB, VIS, and near-IR arms, give spectral resolutions of 
5100, 8800, and 5600 respectively, while for the 1.3/1.2/1.2 arcsec 
slit widths, the resolutions are 4000, 6700, and 4300. For all the 



observations, the UVB CCD was read out binned by a factor of two 
in the dispersion direction. 

The instrument has proven to be useful for the study of gravita- 
tional lensed galaxies due to its complete wavelength coverage and 
high sensitivity throughout the spectral range. Physical properties 
of lensed galaxies can be derived from short exposures as demon- 
strated from comm i ssioning data (Dessauges-Zavadsk y et alj2010t 
IPettini et alj l20ld : IChristensen et alj l20ld) . A wealth of infor- 
mation of the physical properties of galaxies can be gained 
from the full medium-resolution UV- to near-IR spectral cov- 
erage for both lensed galaxie s such as the 8 o'clock arc 
dDessauges-Zavadskv e t al. 20 10 . 1201 1|) or even fainter gamma- 
ray burst (GRB) host galaxies dKriihler et al]|2012l) . With the large 
wavelength coverage, X-shooter is also ideal to determine faint 
galax y redshifts, which is neces sary to construct the lens models 
(e.g. lGrillo&Christensenll201lh . Due to the simultaneous cover- 
age, there are no random offsets in the placement of the slit be- 
tween the three arms, which could be a problem for optical spec- 
tra and near-IR spectra obtained with different instruments. Such 
offsets could introduce differences in the flux levels between the 
UV and near-IR galaxy spectra which are difficult to correct for 
unless the continuum emission is well detected and can be com- 
pared with broad-band photometry. In addition, systematic errors 
in wavelength calibrations could be interpreted as kinematic effects 
such as outflows. X-shooter is therefore ideal for the combined pur- 
pose of determining redshift for faint, lensed sources, and allow a 
characterisation of the physical nature of the lens galaxies. 

Accurate mass models of the lensing clusters are required 
in order to turn them into well-calibrated gravitational telescopes 
because of the degeneracy between the source redshift and lens 
mass. To determine redshifts of lensed galaxies does not require 
high signal-to-noise ratio spectra, and since the UV- selected high- 
redshift galaxies are expected to be star-forming, and therefore 
have emission lines that are much stronger than the underlying 
continuum, a redshift determination can be obtained even for faint 
(m = 25) galaxies with a short integration time. With fairly short 
exposures, we can pinpoint interesting targets for deeper follow-up 
studies. 



2.1 Target selection 

Some of the best calibrated cluster lenses are the z < 0.3 Abell 
clusters. In this paper, we present spectra of galaxies without pre- 
viously determined spectroscopic redshifts lensed by Abell 1689, 
where we aim to determine spectroscopic redshift for galaxies 
which previously only had photometric redshifts. Partly, we focus 
on the second large-scale dar k matter concentratio n in the North- 
East quadrant of the cluster dLimousin et alj|20071) . such that the 
galaxy scale mass substructure can be better constrained in this re- 
gion. 

We also choose to study galaxies lensed by higher redshift 
clusters. We select targets belonging to the MAssive Cluster Sur- 
vey (MAC S), which are massiv e X-ray luminous clusters at 0.3 < 
2 < 0.7 lEbe ling et all 1200 lh . and include observations of the 
southern equivalent to the MACS with declinations < -40 deg (H. 
Ebeling, unpublished). These targets are denoted with the prefix 
SMAC S. Some of the tar gets are also included in the CLASH sur- 
vey dPostman et al.l2012T) . allowing us to investigate in detail galax- 
ies for which high-quality, high-resolution images are available. Al- 
though the 16-broad band filter imaging allows detailed investiga- 
tions of the spectral energy distribution of the galaxies and their 
spatially resolved properties, more accurate information can be ob- 
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tained by spectroscopic follow-up observations. The target galaxies 
are selected from having high total brightness, and from those we 
investigate regions with the highest surface brightness. This im- 
plies that we preferentially select UV bright galaxies that are ex- 
pected to show significant emission lines, because the high surface 
brightness regions are H II regions. Thereby, we can most easily de- 
termine redshifts for the faintest galaxies, and also investigate the 
kinematics and metal abundances in the emission line regions. 



2.2 Spectroscopic observations and data reduction 

The observations were done as a part of the X-shooter consortium 
guaranteed observing time during several different runs in 2010- 
201 1 . A log of the observations can be found in TableQ] The names 
of the arcs targeted for the spectroscopic analysis are hereafter 
shortened, such that the names reflect the cluster name only for 
the MACS clusters, i.e. the arc observed in SMACS J0304.34402 
is named M0304, and the arcs in Abell 1689 are denoted 'A' plus 
by their IDs, e.g. A4.1. In Table [T] and hereafter we adopt the im- 
age IDs from Limousin et al. (2007J) for Abell 1689 i mages, IDs 
for M03 29 in lZitrin et alJd2012al) , IDs for Ml 206 fro m lZitrin etail 
d2012bh . and the ID for M2129 in lZitrin et al] < l201ll) . 

For target acquisition, relative offsets from bright stars were 
used, so the location of the target on the slit depends on the accu- 
racy of the known star coordinates relative to the targets. The slit 
position angle was not preferentially chosen along the long axis of 
the extended arcs. In most cases, we chose a perpendicular angle 
such that nodding the target along the slit would improve the resid- 
uals from sky emission line subtraction. Further, the angle was cho- 
sen to avoid nearby objects. The orientations of the slits are over- 
layed on the panels in Fig. [T] A different procedure was used dur- 
ing the first run, when the Abell 1689 cluster arcs were observed, 
where we opted for longer time on target in stare-mode and did not 
move the target along the slit. In addition, two arcs from different 
sources were inside the slit. Whereas this procedure did allow us to 
determine the arc redshifts, the quality of the near-IR spectra after 
the data reduction was dominated by errors in sky subtraction, in 
particular near sky emission lines. No separate sky exposures were 
obtained. 

The data was r e duced with the ESQ p ipeline version 1.3.7 
dGoldoni et all 120061 : iModigliani etail 120101) using the physical 
model reduction chain teristow et al.ll2008h in both nodding and 
staring mode. The data reduction works in a similar way for the 
three arms. First, the bias level is subtracted from the UVB and VIS 
data (dark current in the case of the near-IR data). The position of 
the orders on the detectors were traced and a two-dimensional flat 
field was created. Then a two-dimensional wavelength solution was 
determined from calibration frames taken the day after the observa- 
tions. Cosmic rays wer e rejected from the science frames using the 
LA Cosmic procedure ( van Dokkum 2001) within the pipeline. Fi- 
nally, the echelle orders were extracted, rectified and merged onto 
a two-dimensional spectrum using the errors as weights for the re- 
gion of the overlap between the different orders. 



2.3 Further data processing 

The data from the pipeline needed a few more processing steps 
before a full analysis could be done. To determine spatial off- 
sets along the slit between each exposure, one could rely on FITS 
header information, or measure offsets directly from the data. We 
chose the latter, but because the continuum emission was very faint, 



this was only possible after binning the 2D data as described in 
section 12.3.1 1 Individually reduced two-dimensional spectra were 
combined using their errors as weights. 

One-dimensional spectra were extracted by co-adding rows 
of data. For each object, the apertures were chosen visually from 
the binned 2D spectra with an aperture size matching for the three 
arms. Typically, an aperture size of 1-2" was used. Source 8.1 in 
Abell 1689, where the X-shooter slit was oriented along the arc, 
was noticeably more extended in the 2D spectrum, and an aperture 
of 4" was used. 

Flux calibration was performed with standard procedures in 
IRAK Observations of various spectrophotometric standard stars 
(GD71, LTT3218, and EG 274) with a 5" wide slit taken at the 
night of the observations were used to calibrate the transmission in 
each of the arms. Reference fl uxes for these stars cov ering the UV 
to the near-IR were taken from Bohlin et al. ( 2001) and Vern et et"al] 
(2008). A correc tion for an avera ge atmospheric extinction mea- 
sured at Paranal (Patat et al]201ll) was applied to the data. M1206, 
M1311, M2031, and M2129 were observed in non-photometric 
conditions. Since the fluxes in the X-shooter spectrum were cor- 
rected to absolute values by the magnitudes measured in HST im- 
ages, flux calibration errors would be partly mitigated as the line 
fluxes were corrected for slit losses. Even though the slit orientation 
was not equal to the parallactic angle, differential slit losses should 
be minor effect, since X-shooter has an atmospheric dispersion cor- 
rector in the UVB and VIS arms. The near-IR arm does not have 
a corrector for atmospheric dispersion, but since the observations 
were done at airmasses below 1.6 the effect was also minor. The 
X-shooter pipeline gives as an output the wavelength solution mea- 
sured in air, which we correct to vacuum heliocentric frame before 
deriving redshifts. Finally, the sp ectra were corrected for Galactic 
reddening using the dust maps of Schle gel et al.l |l998). 

The extracted and flux-calibrated spectra show small offsets 
in the continuum flux between the different arms. In particular, the 
VIS to the near-IR differ by 10-20% in some cases, but given the 
low signal-to-noise ratios of the spectra, the flux levels were not ad- 
justed. The continuum emission follows the expectations over the 
entire wavelength range as illustrated in Fig. [2] In addition to these 
small offsets, the flux calibration showed increased errors at the 
edges of the spectral ranges in each arm, where the total transmis- 
sion is very small due to the dichroic crossover between the arms. 
The spectra in these affected regions are excluded in Fig. [2] 

Emission lines were identified by visual inspection, and line 
fluxes were measured with ngaussf it within IRAF, assuming 
that the line profiles have a Gaussian shape. In some cases of 
very strong lines, the lines are clearly not Gaussian, but can be 
fit well with a combination of two Gaus sian functions (see also 
Pettini et al. 2010; Christensen et al. l2010h . In those cases the listed 
fluxes correspond to the sum of the two components. Uncertain- 
ties of the line width, height, and continuum level are propagated 
for the derivation of emission line flux uncertainties. Some of the 
faintest lines have a more uncertain determination of the peak emis- 
sion wavelength, which implies an apparent shift in redshift com- 
pared to the brightest lines. The shifts can be up to ~100 km s~ , 
but are unphysical if we assume that the lines arise in the same 
medium. We determine source redshifts from the mean weighted by 
the line flux, which overcomes the problem of the small variation 
of redshift. Balmer emission lines are affected by underlying stellar 
absorption. In order to correct for this, we subtracted a stellar pop- 
ulation model described in Sect. 14. 1 1 from the galaxy spectrum to 
get a pure emission line spectrum before fitting the emission lines. 

Removal of telluric absorption lines for the low signal-to- 
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Figure 1. Thumbnail HST images of the clusters around the targeted sources. The images were taken with ACS, WFPC2 or WFC3 with either the F606W, 
F8 14W or Fl 10W filters, and obtained through the MAST or HLA data archive jhttp : //hla ■ stsci ■ edu/fr . The images have varying sizes between 40 
and 90 arcsec on a side with the orientation North up and East left, and with the X-shooter 11" slit position indicated. 



Table 1. Observing log 



Cluster 


Arc ID 


RA (J2000) 


Dec (J2000) 


Slit angle 


Slit widths 


UVB/ VIS/ NIR exposures 


Date 










(deg) 


(") 


(s) 






Abell 1689 
Abell 1689 


4.1 
8.1 


13:11:32.17 
13:11:32.30 


-01:20:57.3 
-01:20:50.8 


-19.4 


1.3/1.2/1.2 


2x3000/2x3000/6x520 


22-03- 


-2010 


Abell 1689 
Abell 1689 


9.1 
31.1 


13:11:30.30 
13:11:30.42 


-01:19:48.33 
-01:19:51.5 


18.8 


1.3/1.2/1.2 


2x3000/2x3000/6x520 


22-03- 


-2010 


SMACS J0304.3^1402 


1.1 


03:04:20.29 


^14:02:27.8 


90 


1.0/0.9/0.9 


4x900/4x900/4x900 


30-08- 


-2011 


MACSJ0329.6-0211 


1.2 


03:29:40.18 


-02:11:45.8 


112 


1.0/0.9/0.9 


4x1200/4x1200/4x1200 


30-08- 


-2011 


SMACS J0359.2-7205 


1.1 


03:59:13.46 


-72:05:14.6 


49 


1.0/0.9/0.9 


3x900/3x900/3x900 


30-08- 


-2011 


MACS J0416. 1-2403 


1.1 


04:16:09.85 


-24:03:42.5 


-160 


1.0/0.9/0.9 


8x800/8x770/8x900 


13-01- 


-2011 


MACS Jl 115.8+0129 


1.1 


11:15:51.40 


01:29:36.6 


-30 


1.0/0.9/0.9 


4x800/4x770/4x900 


12-01- 


-2011 


MACSJ1206.2-0847 


2.1/3.1 


12:06:14.47 


-08:48:33.3 


40 


1.3/1.2/1.2 


4x1070/4x1100/4x1200 


03-04- 


-2011 


MACS Jl 3 11. 0-03 10 


1.1 


13:11:01.37 


-03:10:50.3 


90 


1.0/0.9/0.9 


4x800/4x770/4x900 


03-04- 


-2011 


SMACS J203 1.8-4036 


1.1 


20:31:52.89 


^10:37:32.6 





1.0/0.9/0.9 


4x800/4x770/4x900 


11-10- 


-2010 


MACS J2129 .4-0741 


1.5 


21:29:26.07 


-07:41:41.6 





1.0/0.9/0.9 


4x800/4x770/4x900 


11-10- 


-2010 



noise ratio data increases significantly the noise in the data. Instead 
of correcting the full one-dimensional spectra, which would intro- 
duce additional errors in the spectra, we chose to correct only the 
emission line fluxes individually when necessary. To correct for this 
absorption, we used data of bright, hot stars, observed immediately 
before or after the science exposure and at the same airmass. The 
instrument setup was the same, and the data were reduced in the 
same way as the science frames. A one-dimensional spectrum of 
the star was extracted, normalised, and the ID science spectrum di- 
vided by this normalised spectrum. Regions where the atmospheric 
transmission was less than 5% was set to zero, and when fitting 
emission lines, these regions were excluded. 

2.3. 1 Optimal 2D spectral binning 

Observations in the near-IR from the ground are heavily affected by 
the strong sky emission lines which vary with time. However, blue- 
wards of 2.2 micron, the continuum emission from the sky between 
the sky lines presents a minor contribution. This is the purpose of 
the design of X-shooter; at a sufficiently high spectral resolution, 
one can eliminate the regions of sky lines, and co-add the remaining 
part of the spectrum to obtain a much deeper detection limit than 
would otherwise have been possible as long as the detector read 
noise does not dominate over the background sky emission. 

With this basic principle, we can use the X-shooter data to 
detect very faint objects. In practise, we take the 2D spectra and 
the associated 2D error spectra and bin the data by any required 
amount using a weighted mean algorithm. This binning effectively 
smoothes out any spectral characteristics, so depending on the type 
of investigation, we choose a binning factor between 1 A/pixel (for 
investigating absorption lines and very weak emission lines), and 
~ 200 A/pixel (for investigating the continuum emission from the 
faintest targets). The latter flux calibrated one-dimensional spectra 
for the arcs are presented in Fig. [2] 

2.4 High resolution images 

Since the targets are lensed by well-studied massive galaxy clus- 
ters, most of these have extensive high spatial resolution images al- 
ready available in the HST archive. Both deep and snapshot images 
with the WFPC2, ACS or WFC3 instruments have been taken. We 
obtained these images f rom the HST archive, and used SExtractor 
( iBertin & Arnoutsll 19961) to derive the photometry of the sources as 
listed in Table [2] Up to date zeropoints were used, and the fluxes 
were corrected to total magnitude measured in the detection band 
(/ band). For the most extended arcs SExtractor was first used to 
subtract the extended emission from brighter galaxies in the field 
to create a residual image free of emission from the lens galaxies. 
Then we defined a polygon around the extended arc in the deep- 
est image and derived the magnitude in the other filters using the 
same polygon as reference. The photometric uncertainties were de- 
termined from from photon noise measured in the original images 
before drizzling. In the case of Abell 1689 ID 9.1, the photometry 
was defined for a point source as the 3a limiting magnitude in the 
measured aperture. 

2.4.1 Slit losses 

Besides the loss of flux in slit spectra since some emission fall out- 
side the slit due to the seeing effects and inaccuracy of the center- 
ing on the slit, the lensed images also experience slit losses due to 
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Table 3. Lens image slit loss- and image magnification factors, fi 



Lens image 


ID 


/slit 


ft 


Abell 1689 


4.1 


2.7±0.1 


23.0±3.4 


Abell 1689 


8.1 


2.9±0.7 


38.9±4.0 


Abell 1689 


9.1 


1.0±0.1 


12.7±2.8 


Abell 1689 


31.1 


1.7±o.l 


26.6±3.1 


SMACS0304 


1.1 


3.0±0.1 


42.0±8.0 


SMACS0359 


1.1 


3.0±0.1 


18.0±6.0 


MACS0416 


1.1 


2.6±0.3 


4.6±0.2 


MACS 11 15 


1.1 


5.0±1.7 


16.3±1.4 


MACS 1206 


2.1 


5.1±0.2 


9.2±0.6 


MACS1311 


1.1 


12.7±5.3 


28.0±12.0 


SMACS203 1 


1.1 


4.2±0.3 


15.8±7.0 


MACS2129 


1.5 


2.1±0.1 


9.4±2.0 



their extended and arc-like nature. Some of the observed images 
appear as compact sources even though they are stretched by the 
lensing effects, while other images are very extended low surface 
brightness arcs. Therefore we determine in each case the fraction 
of light falling within the slit relative to the total flux in each band. 
The flux-calibrated spectrum is multiplied with the respective filter 
transmission function to derive the magnitude and the slit losses are 
calculated for each filter. The final, wavelength-independent correc- 
tion factor to be applied to the spectra, is the weighted average of 
the factors in all measured bands, / S ut, as listed in Table [3] where 
the weight in each band determined from the sum of the square of 
the photometric uncertainties in Tableland the uncertainty derived 
from the associated one-dimensional error spectrum. 

Properties derived from spectra obtained in non-photometric 
conditions are partly recovered from the correction of slit losses, 
since the method we adopt does not discriminate between losses 
of light falling outside the slit or loss of light due to absorption 
by clouds. However, since we only apply a constant offset to the 
flux calibration, a wavelength-dependent error caused by the non- 
photometric conditions could be introduced. Also, galaxies have 
colour gradients, but we do not include corrections based on these 
gradients in this investigation. 



2.5 Lens models 

High spatial resolution images are crucial information for the anal- 
ysis of gravitational lensing. The relative projected positions of the 
lens and source galaxies along with spectroscopic (or photometric) 
redshifts are necessary ingredients to determine the total mass dis- 
tribution of the lenses. The lens models allow us to reconstruct the 
original images of the galaxies in the source plane before lensing, to 
derive their radius, which is essential to derive dynamical masses. 

In this paper, we use mass m odels constructed using the 
Lenstool software dJulloetal.ll2007Fl based on the multiple im- 
ages observed in the HST im ages. The detailed Len stool model of 
A 1689 has been published bv lLimousin et alj d2007t) and the mod- 
els for the MACS clusters will be presented elsewhere (Richard 
et al. in preparation). We summarise here the main ingredients 
of the models. Our starting point is the set of multiple systems 
identified in the HST images (between 1 an d 7 systems) which 
complete the ones presented by other groups (Ebelin sTet alj|2009l ; 



1 publicly available at http://lamwws.oamp.fr/lenstool/ 
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Figure 2. Binned source spectra with the best fit spectral templates overplotted by the red lines (Sect. PTTV The grey lines are the error spectra. Regions with 
very noisy data are omitted, both in the fitting and in the display. In particular, wavelengths at 1 fim between the VIS and the near-IR arm are dominated by 
uncertainties in the flux calibration. See the electronic issue of the paper for a colour version of this figure. 
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Figure 2. - continued. The continuum emission in the near-IR for the Abell 1689 arcs, is considerably more noisy compared to the other spectra due to 
increased sky subtraction errors. 



IZitrin et alj|201 ll i lZitrin et alj|2012bll3) . The majority of these sys- 
tems have been confirmed with a measured spectroscopic redshift, 
including the ones obtained from X-shooter and presented in this 
paper. We use the positions and redshifts of these systems to con- 
strain a multi-component mass distribution described as a paramet- 
ric model , including both clu ster-scale and galaxy-scale compo- 
nents (see lRichard et alj|201(t for more details). The Monte-Carlo 
Markov Chain sampler created by Lenstool provides us with a fam- 
ily of best models fitting the constraints, allowing us to derive for 
each parameter its 1 a error bar. The estimate on the magnifica- 
tion factors and their corresponding errors are provided for each 
target in Table [3] In the case of SMACS0304, the uncertainty in 
the cluster redshift (estimated between 0.3 and 0.5) is taken into 
account when deriving the magnification factor of the target 1 . 1 in 



this cluster. With the information of the total flux from the galaxies 
including the effect of magnification and slit losses, we can derive 
the absolute properties of the galaxies. 



3 SOURCE SPECTRA 

In this section we present notes for each of the observed arcs. Ta- 
ble[4]summarises the measured redshifts for the arcs from emission 
and absorption lines, as well as redshifts derived from breaks in 
highly binned spectra, when no clear emission or absorption lines 
could be distinguished. Individual tables in this section list the de- 
tected emission lines and their fluxes uncorrected for magnification 
or slit losses. Typically, only few absorption lines are found by vi- 
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Table 2. Total AB magnitudes of the arcs. 



Cluster 


Arc ID 


F475W 


F606W 


F625W 


F775W 


F814W 


F850LP 


Camera 


A 1689 


4.1 


24.29±0.11 


— 


24.06±0.10 


23.78±0.08 


- 


23.42±0.13 


ACS 


A1689 


8.1 


23.29±0.06 




22.69±0.04 


22.45±0.04 




22.40±0.08 


ACS 


A1689 


9.1 


>28.2 


- 


>27.9 


25.75±0.09 


- 


25.89±0.15 


ACS 


A1689 


31.1 


24.02±0.07 


- 


24.27±0.12 


24.08±0.11 


- 


24.74±0.29 


ACS 


M0304 


1.1 


- 


20.19±0.06 










ACS 


M0359 


1.1 




21.58±0.05 


- 


- 




- 


WFPC2 


JVJ.U £ t 1 


1 1 
1 . 1 




71 T7_|_n HQ 






7i_i_n 1 o 

Zj. 1 j±\)AZ. 




WrrLz 


M1115 


1.1 


23.00±0.41 


23.80±0.10 


22.94±0.37 


23.69±0.40 


23.40±0.6 




ACS 


Ml 206 


2.1 


24.91 ±0.05 


23.16±0.05 


24.24±0.05 


23.12±0.02 


22.99±0.01 


23.19±0.03 


ACS 


M1311 


1.1 




22.58±0.13 










WFPC2 


M2031 


1.1 




21.19±0.10 










ACS 


M2129 


1.5 








23.48±0.14 


23.04±0.09 


22.69±0.15 


ACS 






F105W 


F110W 


F125W 


F140W 


F160W 




Camera 


M0329 


1.1/1.2 




25.0±0.3 










WFC3 


Ml 115 


1.1 




22.94±0.33 


22.86±0.20 


22.54±0.18 


22.20±0.22 




WFC3 


Ml 206 


2.1 


23.13±0.02 


23.13±0.02 


23.21±0.02 


23.12±0.02 


23.10±0.02 




WFC3 


M2129 


1.5 


21.50±0.05 


21.05±0.05 


20.72±0.05 


20.33±0.05 


19.97±0.01 




WFC3 



Table 5. Emission lines from the A 1689 arc ID 4.1 



line 




Aobs 


z 


flux 6 


[Oil] 


3727.09 


8066.41 


1.16427 


11.8±0.2 


[Oil] 


3729.88 


8072.44 


1.16426 


15.6±0.2 


[Ne III] 


3869.84 


8375.75 


1.16436 


3.7±0.5 


H5 


4102.92 


8880.27 


1.16438 


1.8±0.6 


Hi 


4341.69 


9396.41 d 


1.16423 


5.7±1.1 


[Oiii] 


4960.29 


10735.74 c 


1.16434 


15.4±1.9 


[Oiii] 


5008.24 


10840.01 c 


1.16444 


39.8±2.4 


Ho 


6564.63 


14208.28 d 


1.16437 


51.5±4.4 



Notes. 

a Rest-frame vacuum wavelengths. 

6 Emission line flux in units of 10 — 18 erg cm" 

c Close to sky line. 

d Affected by strong telluric absorption lines. 



sual inspection for some of the galaxies, because the S/N level in 
the continuum is low, and the absorption redshift corresponds to 
the average of the detected lines. As a consequence of the low S/N 
levels, the uncertainties of the absorption line redshifts are higher 
than those based on emission lines. 



3.1 Abell 1689 arc ID 4.1 

Two different lensed images (source IDs 4.1 and 8.1) were placed 
withi n the X-shooter slit. The names of the sources are adopted 
from lLimousin et aD d2007l) . Source ID 4. 1 is reported to have z = 
1.1648 dJullo et alj|2010h 7 and is located at the lower edge of the 
slit during the observations, and sky subtraction errors are high in 
the near-IR arm. From the flux-weighted average of the emission 
lines in Table[5] we measure the redshift z = 1.1643 ± 0.0001. 

Among the other strong optical emission lines that fall 
within the spectral range, H/3 is behind a strong sky-line, the 
[Nil] AA6548,6583 and [S II] AA6717,6730 doublets are affected 
by telluric absorption lines so their fluxes are below the detection 
limit. 

In the 2D spectrum we detect the emission from a third ob- 



ject located 2" north of A4. 1. An extended early type galaxy is 
located at RA= 13:11:32.14, DEC = -01:20:55.33 (J2000) in the 
ACS image, but it is only partly covered by the slit. After rebin- 
ning the 2D spectrum a clear break is seen at 4400 A, consistent 
with a Balmer jump at the Abell 1689 redshift. In the spectrum of 
the third object we detect the Ca II AA3934,3639 absorption dou- 
blet at z — 0.170 and a single emission line at 7682 A consistent 
with Ha at z = 0.1702. 



3.2 Abell 1689 arc ID 8.1 (Giant arc) 

Despite being a giant arc, the redshift of imag e ID 8. 1 remains to be 
measured accurately. Broa dhurst et alj j2005l) give the best fit pho- 



tometric redshift z phot = 2.63 ± 0.48 whilelHalkola et alj J200fj) 
find Zphot = 3.10 ± 0.89 and iLimousin et all (2007) determine 
z = 2.30 ± 0.21 based on their lens model. 

The object is located at the upper edge of the slit, and is af- 
fected by sky subtraction residuals, which limit our ability to get 
redshift information from absorption lines. A faint continuum is 
detected in the binned 2D spectrum. With a binning factor of 200 
and 400 pixels in the UVB and VIS arms respectively, any weak 
and narrow emission or absorption lines are effectively smoothed 
out. No emission lines are detected above a 3a level significance in 
the unbinned spectrum. 

The binned UVB and VIS spectrum of ID 8.1 is extracted 
within an 4" aperture and displayed in Fig. [2] The flux-calibrated 
spectrum shows no clear Lyman limit break; however we find in- 
dications of a break at 4470 A. If we in interpret this break as the 
onset of the Lya forest, it implies a redshift of z = 2.68 ± 0.04 by 
fitting stellar population synthesis spectra as described in Sect. 14.11 
The lack of a detected Lyman limit break is due to the high noise 
level at wavelengths shorter than 3500 A. 

3.3 Abell 1689 arc ID 9.1 

Source ID 9.1 is recognised as an F625W drop-out in the ACS im- 
ages, and has a photometric redshift of 5.16 ± 0.74 and an AB 
magnitude of 25.9 in the F775W band teroadhurst et alj|2005l) . 
Halko la et all fcOOrj) find a similar photometric redshift 4.97±0.78, 
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Table 4. Cluster and arc redshifts 



Cluster 


^■cluster 


Arc ID 


Zcra 


d 

break 


■2-abs 


A 1689 


0.181 


4.1 


1.1643±0.0001 a 






A 1689 




8.1 




2.68±0.04 




A1689 




9.1 




5.05±0.02 




A1689 




31.1 


l.sssgio.ooos 6 




1.8348±0.0012 


SMACS 10304.3^1402 


0.3-0.5 


1.1 


1.9634±0.0002 




1.9635±0.0005 


MACS J0329.6-0211 


0.45 


1.1/1.2 








SMACS J0359.2-7205 


0.296 


1.1 


1.0272±0.0005 






MACSJ0416.1-2403 


0.420 


1.1 


1.8960±0.0003 






MACS Jl 115.8+0129 


0.355 


1.1 




0.58±0.02or3.51±0.01 




MACSJ1206.2-O847 


0.439 


2.1/3.1 


3.0385±O.O0Ol c 




3.0372±0.0015 


MACSJ1311.0-O310 


0.494 


1.1 


1.1043±0.00002 






SMACS J2031.8-4036 


0.331 


1.1 


3.5073±0.0002 




3.5061±0.0013 


MACS J2129.4-0741 


0.589 


1.5 


1.3630±0.0004 




1.3617±0.0012 



Notes. 

Cluster redshifts are from lAbell etafl Jl989t) and for Ml 1 15 and M1206 in lEbeling et Zl feOlfJ) . and M2129 in lEbeling et Zl M0359 and M203 1 

cluster redshifts are based on ESO/EM MI and FORS arc hive spectra. The redshift of the M 0304 cluster is measured from photometry. 
Other source spectroscopic redshifts: a |jullo et all fcoiO) report z = 1.1648. b Broa dhurst et al.M2005t) report z = 1.83. c Zitri n et alj f2012bT) report 
z = 3.033, d Redshifts based on spectral modelling. 



while lens models pre dict a lower redshift z mo dci = 2.69 ± 0.27 
dLimousin et alj2007t) . 

Two sources (IDs 9.1 and 31.1) were located within the X- 
shooter slit as illustrated in Fig.[T] The spatial offset between source 
IDs 9.1 and 31.1 projected on the sky is 2'.'5 as measured in an ACS 
image. The A9.1 continuum spectrum is faint, and corresponds to 
an AB magnitude of 26.6 ± 0.4 at ~8000 A, but is detected in 
the binned frames illustrated in Fig. [2] where we identify a break 
at 7381 A. If the break corresponds to the onset of the Lya forest 
it implies z = 5.05 ± 0.02. No Lya emission line nor any other 
significant spectral features are present at these faint flux levels. 



3.4 Abell 1689 arc ID 31.1 

The arc ID 3 1 . 1 is adopted from lLimousin etaL who argued 

that the three images of source ID 12 in iBroadhurst et al . ( 2005) 
are inconsistent with being the same galaxy. The source 31.1 was 
located close to the edge of the slit, where sky subtraction errors 
were higher than the average. 

A wealth of emission lines are detected as listed in Table [6] 
Excluding Lya while weighting the other emission line redshift by 
their measured fluxes , we find z = 1 .8339 ± 0. 0003, consistent 
with t hat reported in IBroad hurst et al. I d2005h and iLimousin et al.l 
(2007). The Lya line profile is double peaked with a brighter 
red than blue component and has an equivalent width of 40 
A. We detect several uncommon emission lines for this source. 
In particular we draw attention to the Olll] AA1661,1666 and 
Civ AA1548,1550 doublets, which indicate the presence of a 
strong ionisin g source. These line s were also detected in a nother 
lensed galaxy dFosburv et al]|2003l ; Ivmar-Martin et alj|2004l) . The 
temperature - sensitive [O III] A436 3 line is rarely detected in galax- 
ies at z > 1. 1 Yuan & KewlevH2009T) find this line in another lensed 
galaxy at z = 1.7, also behind the Abell 1689 cluster. In a separate 
paper we presen t a more detailed analys is of the emission from this 
unusual galaxy dChristensen et al.ll2012h . 



Table 6. Emission lines from the A 1689 arc ID 31.1 



line 




^obs 


z 


flux 6 


Lya 


1215.67 


3445.77 


1.8345 


72.7=1=0.7 


CIV 


1548.20 


4386.81 


1.83349 


4.3±0.4 


CIV 


1550.77 


4394.22 


1.83357 


4.7±0.3 


Olll] 


1660.81 


4705.82 


1.83345 


3.2±0.4 


Om] 


1666.15 


4720.63 


1.83326 


7.4±0.3 


[Cm] 


1906.68 


5402.29 


1.83335 


9.1±0.4 


Cm] 


1908.73 


5408.19 


1.83340 


4.7±0.4 


[On] 


3727.09 


10563.06 


1.83413 


9.2±0.8 


[On] 


3729.88 


10569.40 


1.83371 


11.7±0.9 


[Ne III] 


3869.84 


10965.28 


1.83352 


6.4±2.5 


H7 


4341.69 


12301.90 


1.83344 


17.0=1=2.1 


[Om] 


4364.44 


12366.36 


1.83344 


6.2=1=1.8 


H/3 


4862.70 


13777.99 d 


1.83340 


42.7±4.9 


[Om] 


4960.29 


14057.55 d 


1.83402 


63.0±10.6 


[Om] 


5008.24 


14194.09 c ' d 


1.83415 


211.5=1=5.7 



Notes. 

a Rest-frame vacuum wavelengths. 

b Emission line flux in units of 10~ 18 erg cm -2 s . 

c Close to sky line. 

d Affected by strong telluric absorption lines. 



3.5 SMACS J0304.3-4402 arc ID 1.1 

The morphology of the lensed galaxy is complex, and likely con- 
sists of several components that are in the process of a merger. In 
each of the lensed image counterparts in the HST image, five indi- 
vidual regions are visible. The X-shooter spectrum of the source 
includes three of these regions. Due to blending by the seeing, 
however, only two regions are spatially distinguishable in the two- 
dimensional spectrum. Many strong rest-frame optical emission 
lines are detected from this galaxy at z = 1.9634 ± 0.0002 as 
listed in Table [7] including again the O III] AA1661,1666 doublet. 
The emission lines from the two regions have a velocity offsets 
of ~100 km s _1 , while the UV continuum appears similar in 
shape. The emission line fluxes in Table [7] are the total flux in 
the X-shooter slit summed for the two regions. Even though the 
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Table 7. Emission lines from the M0304 arc ID 1.1 



Table 8. Emission lines from the M0359 arc ID 1.1 



line 




Aobs 


2 


a b 

flux 


O III] 


1660.81 


4921.93 


1.96357 


7.1±1.6 


in] 


1666.15 


4937.35 


1.96333 


11.4±1.9 


[On] 


2470.22 


7322.39 


1.96337 


8.4±2.1 


[On] 


3727.09 


1 1044.78 


1.96331 


437.8±3.4 


[On] 


3729.88 


11053.09 


1.96332 


587.2±3.9 


H9 


3836.49 


11370.55 c 


1.96372 


29.1 ±4.2 


[Nelll] 


3869.84 


11468.93 c 


1.96360 


149.9±3.8 


He I + H8 


3890.17 


11527.64 c 


1.96320 


95.3±3.2 


[Nelll] 


3968.53 


11761.36 


1.96359 


64.3±3.4 


H7 


3971.20 


1 1767.48 


1.96313 


51.8±1.7 


Hi 


4102.92 


12158.43 


1.96329 


128.5±2.7 


Hi 


4341.69 


12866.19 


1.96333 


205.9±1.3 


H/3 


4862.70 


14410.72 c 


1.96345 


500.0±0.9 


[OIII] 


4960.29 


14700.41 c 


1.96355 


669.5±1.9 


[OlII] 


5008.24 


14841. 33 c 


1.96331 


2357±3.5 


He I 


5877.28 


17417.44 


1.96345 


47.0±5.2 


Ha 


6564.63 


19453.85 c 


1.96336 


1802.1 ±0.6 


[Nil] 


6585.42 


19513.96 c 


1.96313 


58.1±2.7 


[Sn] 


6718.29 


19910.44 


1.96355 


156.8±1.1 


[Sll] 


6732.67 


19951.76 


1.96335 


102.7±1.7 



Notes. 

a Rest-frame vacuum wavelengths. 
6 Emission line flux in units of 10~ 
c Affected by telluric absorption lines 



18 erg cm 2 s 



Balmer lines are very bright from this source, Lya emission is ab- 
sent, and its (UV) spectrum appears similar to the cB58 spectrum 
I Pettini et"al] 200C ). A m ore detailed investigation is presented in 
IChristensen et all J2012h . 



line 




Aobs 


z 


flux 6 


[Oil] 


3727.09 


7558.09 


1.02788 


97.7±3.2 


[Oil] 


3729.88 


7560.08 


1.02690 


137.9±3.4 


[Ne III] 


3968.53 


8048.88 


1.02818 


5.0±1.7 


H5 


4101.74 


8317.41 


1.02778 


11.0±2.5 


H~, 


4341.69 


8800.95 


1 09708 


37 1±1 6 


H/3 


4862.70 


9856.80 


1.02702 


95.4±1.4 


[OlII] 


5008.24 


10152.15 


1.02709 


87.9±4.0 


[Nil] 


6549.91 


13276.12 c 


1.02692 


23.8±3.8 


Hq 


6564.63 


13306.85 c 


1.02705 


369.1 ±0.9 


[Nil] 


6585.42 


13348.02 c 


1.02690 


55.4±4.4 


[Sill] 


9533.20 


19324.30 c 


1.02705 


56.4±4.0 


Notes. 










a Rest-frame vacuum 


wavelengths. 






b Emission line flux in units of 10 18 


erg cm~ 2 


s- 1 . 


c Affected by telluric absorption lines 






Table 9. Emission lines from the M0416 arc ID 1.1 


line 




•^obs 


z 


flux 6 


[On] 


3727.09 


10793.11 


1.89585 


17.2±1.4 


[On] 


3729.88 


10801.18 


1.89585 


24.5±1.4 


H/3 


4862.70 


14081.70 c 


1.89586 


27.4±2.1 


[Om] 


4960.29 


14367.79 c 


1.89656 


27.2±2.8 


[Om] 


5008.24 


14503.75 c 


1.89598 


100.2±4.2 



Notes. 

a Rest-frame vacuum wavelengths. 

b Emission line flux in units of 10 — 18 erg cm -2 s~ 

c Affected by strong telluric absorption lines. 



3.6 MACS J0329.6-0211 arc ID 1.1/1.2 

The two mirrored lensed images with ID s 1.1. and 1.2, whic h have 
a photometric redshift of z — 6.18^0 05 dZitrin et"aTll2012ah . were 
targeted by the observations, which were carried out by both nod- 
ding and offsetting along the slit. The galaxy has a high surface 
brightness region at the edges pointing towards each other, and a 
more elongated lower surface brightness region in the extended 
arcs. The slit orientation covered only the high surface brightness 
regions, where the brightest image I D 1.2 has a magnit ude of 25.0 
AB mag in the F110W image (see IZitrin et alj|2012al for photo- 
metric data of all lensed images). 

We detect neither emission lines nor the continuum emission 
even in a highly binned 2D spectrum. With the faintness of the 
galaxy, we do not expect to detect the continuum emission even 
with a very large binning factor, but aim instead for a detection 
of a single emission line. We would be able to detect Lya in the 
VIS arm to a 3a significance level of 5 x 10 -18 erg cm -2 s _1 
assuming a typical line width of 50 km s _1 . We cannot exclude 
that an emission line happens to fall directly at the wavelength of 
a strong sky emission line given the uncertainty of the photometric 
redshift. In that case, the 3a detection limit could be much higher: 
~2x 10~ 16 erg cm" 2 s" 1 . 

3.7 SMACS J0359.2-7205 arc ID 1.1 

Emission lines from this arc are listed in Table[8] No unusual strong 
UV emission lines are detected even though its low redshift (z = 



1.0272 ± 0.0005) and the wavelength coverage would allow for 
detection of the [C III],C III] doublet. [Om]A4959 lies in a very 
noisy region where the transmission of the NIR arm is low and is 
not detected. 

3.8 MACS J0416.1-2403 arc ID 1.1 

Emission lines detected from this arc at z = 1.48960 ± 0.0003 
are listed in Table |9] H/3 and the [O in] AA4959,5007 doublet are 
strongly affected by telluric absorption lines, and the correction for 
the atmospheric absorption implies a correction factor between 3 
and 10. The flux uncertainty only reflects the emission line fitting, 
since the error for the correction for telluric lines is not propagated. 
Again, no strong UV emission lines are detected, because even the 
strongest rest-frame optical lines are relatively faint. 

3.9 MACS J1115.8+0129 arc ID 1.1 

The image of the source is a very low surface brightness arc. The 
X-shooter slit was oriented across the arc so the slit losses were 
substantial. 

We do not detect any emission lines in either of the arms. The 
continuum emission is the faintest among our sample, and is only 
recognisable after binning the spectrum by a factor >20 along the 
dispersion. The binned extracted spectrum in Fig. [2] is extracted 
with an aperture of 1"5. The NIR spectrum is also very faint, but the 
continuum emission is detected in the binned spectrum although it 
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Figure 3. Upper panels: Two-dimensional X-shooter spectrum around the 
single emission line detected for the Ml 206 source, and the corresponding 
region around the source in an 7/S77ACS/F606W image. The two source 
IDs 2.1 and 3.1 are separated by 2'.'3 and are seen in both panels. Source 
2. 1 to the North-East in Fig. [T] (upper part in this image) has a compact 
morphology, and source 3.1 to the South- Western one is more extended. 
The one-dimensional spectra in the middle and lower panels are extracted 
for each of the sources and are corrected for telluric absorption lines. 

is too faint to be detected in the K band. Between the UVB and VIS 
spectra, there appears to be a break, which we can interpret either as 
the Balmer decrement at z ~ 0.5, or a higher redshift for the onset 
of the Lya forest at z ~ 3.5. The exact wavelength of the break 
is sensitive to the binning of the data; a smaller binning suggests 
a higher break wavelength, and thus a higher redshift. Leavi ng the 
redshift as a free parameter in HyperZ (Bolzonella et al. 200(J), we 
determine the galaxy redshift z = 0.58 ± 0.02 or z = 3.51 ± 0.01. 
The physical properties of the galaxy derived in Sect. [4] depend 
on the assumed redshift, and we suggest that the lowest redshift is 
valid one based on the lensing model. 

3.10 MACS J1206.2-0847 arc IDs 2.1 and 3.1 

The conditions during the observations were poor with chang- 
ing seeing and thin cirrus clouds. The slit orientation was chosen 
su ch that the two com ponents, denoted by source IDs 2.1 and 3.1 
in Zit rin et ID {2012b), which are separated by 2" 3, were placed 
within the slit. The two sources have the same redshift jZitrin et aU 
l2012bl) . 

We only detect a single, but remarkably bright emission line 
at 20220 A, as illustrated in Fig. [3] and a break in the contin- 
uum around 3700 A. The line is identified as [O III] A5007 at 
z — 3.038. At this redshift, all the other strong rest-frame opti- 
cal emission lines ([O II] AA3727,3730, H/3, and [O III] A4959) fall 
at the wavelengths of sky emission lines or are affected by telluric 
absorption. At weaker significance, we detect absorption lines at 
2 = 3.0372 ±0.0015 in the binned spectra illustrated in Fig.|4] The 
marked lines are the strongest UV lines d etected in other hi gh red- 
shift galaxies, such as in the cB58 galaxy JPettini et alj20 02). It ap- 
pears that Lya is seen in emission at the red wing of a damped Lya 
absorption trough, however, the emission line is less significant, 
and cannot be discerned in the unbinned one- or two-dimensional 
spectra. The inferred redshift is con sistent with that me asured with 
low-resolution optical spectra only jZitrin et alj2012bl) . 

The emission line profile of the compact north-east image ID 
2.1 appears to be double-peaked, while the source ID 3.1 only 
shows one as illustrated in Fig. [3] The separation of the two emis- 
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Figure 4. Section of the UVB and VIS spectra of the M1206 source ID 2. 1 
binned by 10 and 20, respectively, to a dispersion of 2.5 A pixel -1 . Typical 
strong UV absorption lines seen in high redshift galaxies are marked at the 
redshift z = 3.0385. 



Table 10. Emission lines from the Ml 206 arc. Source IDs are adopted from 
lzitrinetall2012bl) . 



line A r a ^obs z flux 6 
Source ID 2.1: 

[O III] 5008.24 20217.87 3.03696 83.4±1.4 

[Om] 5008.24 20225.77 3.03854 101.9±1.8 

Source ID 3.1: [O in] 5008.24 20226.02 3.03859 119.7±1.2 



Notes. 

a Rest-frame vacuum wavelengths. 

b Emission line flux in units of 10~ 18 erg cm~ 2 s — 1 . 

sion components in the 2.1 source is 7.9±0.3 A corresponding to a 
velocity separation of 1 17±5 km s . 

The emission line fluxes in Table 1 101 include the uncertainties 
from fitting Gaussian profiles to the observations, but they do not 
include the uncertainty due to the non-photometric conditions dur- 
ing the observations. The line fluxes are derived after correction for 
telluric lines, which absorb between 40 and 60% of the emission. 



3.11 MACS J1311.0-0310 arc ID 1.1 

From the emission lines listed in Table [TT] we derive a weighted 
redshift average of z = 1.10433 ± 0.00002. At this redshift Ha 
and [Nil] fall in the region between the J and H bands where 
telluric absorption is severe. The observing conditions were non- 
photometric, which affect the absolute values of the emission lines. 
We note that the ratios of emission lines is most likely not affected, 
because the transmission derived from the spectrophotometric stan- 
dard star observed during the night appeared similar in shape to that 
of other nights. 

3.12 SMACS J2031.8^1036 arc ID 1.1 

Several emission lines are detected from this arc as listed in Ta- 
ble [12] and the flux-weighted average redshift is z — 3.5073 ± 
0.0002. Most lines are not affected by telluric absorption, but 
H I A397 1 has a significant correction, and the flux in Tablell2|gives 
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Table 11. Emission lines from the M131 1 arc 



line 




Aobs 


z 


flux 6 


[Oil] 


3727.09 


7843.14 


1.10436 


8.0±1.3 


[Oil] 


3729.88 


7848.85 


1.10432 


12.0±0.4 


H/3 


4862.70 


10232.73 


1.10433 


30.6±2.6 


[Om] 


5008.24 


10539.05 


1.10434 


13.9±2.4 



Notes. 

a Rest-frame vacuum wavelengths. 

6 Emission line flux in units of 10 -18 erg cm -2 s -1 . 



Table 13. Emission lines from the M2129 arc 



line 




^obs 


2 


flux 6 


[Oil] 


3727.09 


8805.93 


1.36251 


8.6±1.2 


[Oil] 


3729.88 


8814.78 


1.36315 


7.5±U 


[Om] 


5008.24 


11834.63 


1.36314 


36.5±1.7 



Notes. 

a Rest-frame vacuum wavelengths. 

6 Emission line flux in units of 10 — 18 erg cm -2 s — 1 . 



Table 12. Emission lines from the M203 1 arc 



line 




^obs 


z 


flux 6 


Lyct 


1215.67 


5480.80 


3.50845 


173.7±0.5 


NIV] 


1486.50 


6700.85 


3.50780 


7.6±0.9 


Om] 


1660.81 


7485.15 


3.50693 


2.9±0.6 


Om] 


1666.15 


7509.46 


3.50707 


8.8±0.7 


[Cm] 


1906.68 


8594.04 


3.50733 


12.0±0.5 


Cm] 


1908.73 


8603.54 


3.50747 


8.4±0.9 


[On] 


3727.09 


16800.18 


3.30759 


14.4±1.6 


[On] 


3729.88 


16811.01 


3.50712 


15.5±1.3 


[Ne III] 


3869.84 


17442.96 


3.50741 


14.8±2.1 


Hl(H7) 


3971.20 


17898.76 


3.50714 


6.0±1.0 


H/3 


4862.70 


21917.47 


3.50726 


37.7±1.3 


[Om] 


4960.29 


22357.68 


3.50733 


61.9±0.9 


[Om] 


5008.24 


22573.68 


3.50731 


205.2±0.5 



Notes. 

a Rest-frame vacuum wavelengths. 
6 Emission line flux in units of 10 _ 
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Figure 5. Section of the VIS and near-IR spectra of the M2129 arc ID 1.5 
binned to a dispersion of 5.2 and 12 A pixel -1 , respectively. Detections 
of Call 3934,3969 lines and absorption lines from the Balmer series of 
hydrogen at z = 1.3617 are indicated. Also a faint emission line from the 
[O II] A3727 doublet is visible in this binned version. 



tween 9300 and 9500 A are corrected for, and several Balmer series 
absorption lines are also detected (see Fig. [5}. 



the flux after the correction is applied. Apart from that, telluric ab- 
sorption and sky emission lines affect our ability to detect H7 and 
H<5, which fall in the region between the H and K bands. 

Again in this case, the Lya emission line profile is double- 
peaked, and a wealth of information can be derived from the 
UV emission lines, which is presented in a separate paper 
dChristensen et alj|2012t) . The redshift listed in Table [T2l for Lya 
reflects the redshift of the bluest wavelength of the red peak as iden- 
tified by visual inspection. 

At the low S/N per pixel in the UVB and VIS arms, we can 
only identify several Lyman forest absorption lines, and Si II A1206 
and SlV AA1394,1403 at z = 3.5. A detailed analysis of the ab- 
sorption lines requires a spectrum with a higher S/N measured in 
the continuum. 

3.13 MACS J2129-0741 arc ID 1.5 

The galaxy cluster MACS J2129.4- 0741 at z = 0-58 9 shows a 
galaxy lensed in six distinct images. iZitrin et al. H20Tl|) predict a 
source redshift of z — 1.0 — 1.5 based on their lens model. The 
galaxy is red, has colours similar to the lens galaxies, and its mor- 
phology is consistent with being an early type galaxy. 

For the X-shooter observations, we targeted the isolated, 
southern-most image of the six. Although it is not the brightest 
of the images, source blending and background subtraction is not 
a problem at this location. In the extracted spectra we find weak 
oxygen emission lines listed in Table|13l These lines imply a red- 
shift z = 1.3630 ± 0.0004. At this redshift, absorption lines from 
Ca II AA3934, 3969 are detected after telluric absorption lines be- 



4 PHYSICAL PROPERTIES FROM THE CONTINUUM 
EMISSION 

In this section we derive physical properties based on the contin- 
uum spectra from the sources. 



4.1 Stellar population fit 

To determine the stellar masses of the galaxies, along with other 
parameters, we use spectral energy distribution fitting (SED) meth- 
ods. Conventionally, SED fits make use of sparsely sampled photo- 
metric data from broad or narrow band imaging data, whereas we 
here are dealing with the entire spectral range from the observed 
UV to near-IR. Since the spectra have low S/N ratio, we first bin 
the data by a factor between 100 and 1000, depending on the ini- 
tial quality of the spectra, giving resolutions in the range ~ 20-200. 
The binned spectra and the associated error spectra are shown in 
Fig.H 

We compare the spectra with a large set of model spectra 
from Bruzual & Charlo t H2003h . We extract spectral templates with 
0.2, 0.4 and 1.0 solar metallicities, and star formation histories de- 
scribed as an instantaneous burst or exponentially declining star 
formation rates (SFR) with e-folding time scales of 10, 50, 100, 
300, 500, and 1000 Myr, and exponentially increasing SFRs with 
time scales of 10, 50 and 500 Myr. The initial mass functi on (IMF) 
is eith er a Salpeter (Salpete^ [T955l) . or a Chabrier IMF dChabrieJ 
l2003h with lower and upper mass cutoffs of 0.1 and 100 Mq, re- 
spectively. The best fits are determined by minimising 
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X — /^((/obs.A — /model, a) /fobs, a) , 
A 

where the sum is taken over all photometric points, A in the spectra. 
In practise, there is no significant distinction between the quality of 
the fits between the two IMFs. We therefore choose to only use 
the Chabrier IMF to derive the stellar masses for consistency. If 
a Salpeter IMF is selected instead, the stellar mass is a factor 1.7 
higher due to the flattening of the Chabrier IMF below 1 Mq. 

To determine the best fit spectral template, the flux at each 
(binned) wavelength is converte d to an AB magnitude, and we use 
a modified version of HyperZ (Bolz onella e"tal]l200rjh . which al- 
lows us to include for each object a catalogue which contains up 
to 500 individual photometric data points. From all the photomet- 
ric points between 3000-22000 A, we exclude the wavelengths 
around 1 micron, where the flux calibration causes inaccuracies at 
the wavelengths where the dichroic splits the light between the VIS 
and near-IR arms. We also exclude the regions between the J-H 
and H-K bands which are strongly affected by telluric absorption. 
HyperZ also allows to determine the best fit reddening as an addi- 
tional parameter. He re we use a starburst galaxy extinction curve 
dCalzetti et al. 2000) as input. We do not constrain the reddening 
with that determined from the Balmer decrements where available 
(see Sect. I5.2K because the stars may not experience t he same 
amoun t of reddening as the gas. For starburst galaxies, ICalzettil 
finds that the radiation from stars and ionised gas experi- 
ence different extinctions and that the stellar reddening is on the 
average 44% of that in the gas phase. The flux attenuation blue- 
wards of Lyo is taken into account using the prescriptio n for the 
transmission in the intergalactic medium in Madau) d 19951) . 

For the galaxies with known spectroscopic redshifts, their 
redshifts are fixed in the template fitting, while in two cases 
(MACS 1115 and Abell 1689 ID 8.1), the best fit was also used to 
derive the photometric redshift. We do not constrain the template 
metallicities to the ones best matching the nebular abundances de- 
rived in Sect. 15.41 because these may not be the same, as older stel- 
lar populations, if present, may have lower metallicities reflecting 
the lower abundances of the gas at their time of formation. The out- 
put gives the best fit template which allows us to derive the ages, 
metallicities, and reddening as listed in Tablell4l The best fit spec- 
tra are overplotted on Fig. [2] One of the outputs of HyperZ is a 
scaling factor b allowing the galaxy stellar mass to be derived: 

M, = 6(1 + z)AtyDI/Lq ,bol, 

where Lq^oI is the solar bolometric luminosity (3.826xl0 33 erg 
s _1 ) and Dl is the luminosity distance, where we assume a flat 
cosmology with J7a = 0.72 and Ho — 73 km s _1 Mpc -1 . The de- 
rived stellar mass is magnification dependent and is sensitive to slit 
losses. Tablell4llists the corrected stellar masses including the cor- 
rection factors in Table [3] Uncertainties for the parameters (mass, 
age, reddening) are determined by adding to the initial spectrum 
the error spectrum multiplied by a random number normalised to 
a Gaussian distribution with a — 1. A thousand fits are performed 
for each galaxy spectrum, and the standard deviation of the fitted 
values are taken as the uncertainties for the parameters. 

In all cases the best fit templates are obtained with a model 
of an instantaneous burst of star formation. Such models are ar- 
tificial in the sense that this would mean that there is no current 
on-going star formation, which is not correct because the galaxies 
show strong emission lines. Other templates that produce accept- 
able fits are those with a small value of r — 10 Myr, while in a 
few cases, templates extracted from a constant star formation rate 
is also a valid solution. In a few cases, an exponentially increasing 



SFR also produces an acceptable fit. With these alternative fits the 
stellar masses remain the same as those listed in Tablell4lto within 
±10-20% uncertainty. 

As usually the case for SED fitting, there is a strong degen- 
eracy between the age of the stellar population and the reddening. 
For any given spectrum, a similar good fit can be obtained with a 
younger population combined with more reddening of the template. 
The galaxy ages as obtained by the best fits are generally small 
(< 100 Myr), and the reddenings are small as well, such that the 
allowed parameter space of reddening-age values is limited. One 
concern is that, if the ages are underestimated by the SED fit, the 
stellar masses will be as well. It is common in the literature to only 
include templates in the SED fits that are older than the typical dy- 
namical timescale f or galaxies, i.e. older than a few 10 Myr (e.g. 
IWuvts et al.ll2012bl) . For the lensed galaxies in this study, the re- 
sulting fits when forcing the age to be older than 20 Myr, the fits 
deteriorate significantly bluewards of the Lya wavelengths. Simi- 
larly, if we fit the spectra with templates of a constant star formation 
rate and a lower limit for the age of a few times 10 Myr, the fits at 
UV wavelengths deteriorate, while the rest of the spectrum can be 
fitted well with an older template. We therefore choose to not con- 
strain the possible template ages in the fitting. Another issue is that 
an underlying older stellar population can add to the total stellar 
mass, while not be detected because it does not contribute signifi- 
cantly to the UV radiation. Since the JC-band fluxes of the galaxies 
are small, we argue that a significantly more massive older stellar 
population is not present. To better constrain the exact masses of an 
underlying older stellar populations in the galaxies, deep rest-frame 
near-IR images are needed to fit the SED with a larger wavelength 
range. 

4.2 Absolute magnitudes 

A relevant question to ask given the heterogeneous sample, is what 
their magnitudes are relative to a typical galaxy at the same redshift. 
To calculate this, we use the best fit model spectrum to derive the 
magnitudes at rest-frame 1500 A, and co mpare these to the v alues 
of MXbi which depends on the redshift dArnouts et al.ll2005h . We 
use the best fit templates, rather than the observed spectra because 
of the low S/N ratio in the spectra. The resulting absolute magni- 
tudes and that compared with MXb{ z ) are listed in Table [T5l The 
errors reflect the uncertainties in the magnification factor and slit 
losses. 

For the galaxies with a low UV-flux characteristic for an 
evolved stellar population, neither their UV absolute magnitude nor 
their UV SFRs derived in Sect. 15.51 are relevant for comparing 
with star forming galaxies at a given redshift selected in a uniform 
way from large surveys. When the spectra show evidence for older 
stellar populations it is more reasonable to compare the relative lu- 
minosities at rest-frame visual wavelengths. Table Q3] lists the ab- 
solute magnitudes and luminosities of the galaxies at rest-frame 
6500 A corrected for magnification and slit losses. Again the best 
fit template is used to determine the luminosities. For the Abell 
1689 targets with poor near-IR spectra, and M2031 at z — 3.5, the 
6500 A flux is based on the best fit template spectra at wavelengths 
redwards of the measured spectrum. The redshift dependent char- 
acteristic magnitudes (M£b) defined by the luminosit y function of 
galaxi es in the rest-frame r' band are adopted from Gabas ctTet al.l 
( l2006h . 

Comparing the relative rest-frame UV and optical luminosi- 
ties, we see that most galaxies are UV bright. This is by selection, 
since the observations targeted galaxies with both high total bright- 
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Table 14. Best spectral fit parameters from Sect. |4-1| The masses, SFRs, and magnitudes are corrected for slit losses and lens magnification factors. 



Arc 


Z/Z 


Age 


E(B - y) s tar S 


log M* 


UVSFR 


UV SSFR 






(Myr) 


(mag) 


(M ) 


(M Q yr- 1 ) 


(Gyr- 1 ) 


A4.1 


0.2 


4±3 


0.10±0.02 


6.97±0.08 


0.4±0.1 


40±11 


A8.1 


0.2 


4±3 


0.07±0.08 


7.85±0.15 


1.6±0.7 


23±11 


A9.1 


1.0 


90±50 


0.00±0.02 


9.79±0.13 


0.6±0.2 


0.09±0.04 


A31.1 


0.2 


4±2 


0.17±0.04 


7.70±0.10 


1.0±0.2 


20±6 


M0304 


0.4 


57±6 


0.17±0.02 


10.57±0.08 


11.6±2.3 


0.30±0.02 


M0359 


0.4 


110±18 


0.15±0.03 


9.62±0.13 


1.5±0.5 


0.36±0.05 


M0416 


1.0 


32±14 


0.00±0.08 


8.70±0.16 


0.9±0.2 


1.8±0.8 


M1115 


1.0 


400±260 


0.00±0.06 


7.83±0.22 


0.007±0.003 


0.10±0.06 


Ml 206 


1.0 


20±13 


0.00±0.03 


9.23±0.25 


3.7±0.6 


2.2±1.7 


M1311 


1.0 


140±20 


0.15±0.09 


9.14±0.21 


0.3±0.2 


0.2±0.1 


M2031 


0.2 


4±3 


0.00±0.03 


9.16±0.21 


18.0±8.2 


13±6 


M2129 


1.0 


4700±1500 


0.07±0.06 


11.39±0.12 


0.003±0.001 


10~ s 



ness and high surface brightness regions, which at higher redshifts 
implies a high UV surface brightness. The only exception in our 
sample is the M2129 arc, which was known to be an old evolved 
system based on its red colour and its morphology. 



5 PHYSICAL CONDITIONS FROM EMISSION LINES 

In this section we present properties of the lensed galaxies, which 
have sufficiently bright emission lines for a derivation of physical 
quantities. Diverse characteristics of the galaxies are expected, be- 
cause the galaxies are not selected with a known redshift, nor ap- 
pearance, hence we cannot derive a uniform set of properties for 
all sources with the current data set. All the physical properties 
derived in this section are summarised in Table [T6l When line ra- 
tios are considered, we do not apply any correction for slit losses 
and magnification factors, since we assume that the emission lines 
come from the exact same region. 



5.1 Kinematics 

5.1.1 Line widths 

We derive the mean velocity dispersion (a) from all the emission 
line widths, excluding Lya when detected, and weighted with their 
associated uncertainties. The measured line widths are corrected 
for instrumental resolution by measuring the widths of nearby un- 
blended sky lines, and subtracting the widths in quadrature. The 
results are listed in Table [T6l To determine the dynamical masses 
of the galaxies, their sizes need to be measured after lens model- 
ing and source reconstruction, and will be presented elsewhere (J. 
Richard in prep.). Although not included in Table|16|because only 
few lines are detected, we can still determine the velocity disper- 
sion of the emission lines from M2129 ID 1.5 to be 147± 17 km 
s~\ and 40±3 km s" 1 for M1206 ID 2.1. The small velocity dis- 
persions derived from the emission lines agree with the small stellar 
masses determined from the SED fits. 



5.1.2 Velocity offsets 

High-redshift galaxies show velocity offsets between Lya emis- 
sion and absorption lines, which are generally interpreted as galaxy 
scale outflows caused by feedback from star formation. Detailed 



studies have used the Lya line to measure the emission line red- 
shift in comparison to rest-frame UV absorption lines arising in 
the interstellar medium (IS M) of the LBGs and f ound relative off- 
sets of several 100 km s _1 dShaplev et alj|2003l) , but also fainter 
galaxies identified by GRBs appear to have large velocity shifts 
dMilvang-Jensen et alj|2012h . In a sample of 8 Lya selected galax- 
ies velocity offsets of I45I23 ' fln s_1 between the systemic red- 
shifts from rest-fr ame optical emission lines and the Lya lines have 
been measured ( Hashimoto et al.ll2012l) . These small velocity off- 
sets are similar to the offsets determined from the two Lya emit- 
ting galaxies in this study. Since in the reference frame of the out- 
flowing gas, photons bluewards of the Lya line centre are shifted 
into reson ance, the photons escap e only after having diffused to the 
red wing. IVerhamme et al.l 112006) find that Lya experiences about 
twice the velocity shift of the true outflow velocity in low hydrogen 
column density LBGs. 

Some Lya lines have a prominent blue cutoff in the spectral 
profile while others show both a blue and a much stronger red com- 
ponent. Observed at lower spectral resolution, the intrinsic red- 
shift of the Lya lines is difficult to determine. A better approach 
is to compare the ISM lines with other emission lines than Lya, 
preferably measured with the same instrument to avoid systematic 
offsets. X-shooter spectra of the 8 o'clock arc show that the UV 
absorption lines and optical emission lines have redsh ifts consis- 
tent with each other (Dessauges-Zava dskv et "all l201ll) , but since 
the broad ISM absorption lines span a range of velocities relative 
to the systemic one, a galaxy scale outflow of 120 km s _1 is likely. 
When outflows have been observed in galaxies, their ISM absorp- 
tion lines are blueshifted with respect to the systemic redshift. 

Table [4] lists both redshifts derived from emission and ab- 
sorption lines. The absorption line redshifts are measured by fit- 
ting Gaussian functions and are less accurate because the S/N per 
pixel is low. Only for the sources with the brightest UV continuum 
can we measure reliably the redshifts from absorption lines. Even 
though the continuum is relatively bright, we find no significant 
ISM lines in the M0359 spectrum since at z — 1.0 most strong 
lines are bluewards of the UVB spectrum. Strong low-ionisation 
lines that are relatively easily detected and are not blended include 
Sill A1260, Sill A1304, ClI A1334, and Sill A1526, and among 
the high-ionisation species we detect S IV AA1393,1402 and C IV 
AA1548,1550. LBG spec tra show velocity offsets between high 
and low-ionisation lines ( Shapley et al. 2 0031) . as well as offsets 
between emission and absorption lines indicative of galaxy scale 
outflows. Only in M0304 do we detect several high- and low ioni- 
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Table IS. Rest-frame optical absolute magnitudes and relative luminosities at 1500 and 6500 A. 



Arc 


-Mab,15uo 


L/L*(z,l500A) 


A^AB,6500 


L/L* (2,6500 A) 


A4.1 


-18.39±0.29 


0.2±0.1 


-17.93±0.29 


0.013±0.004 


A8.1 


-20.56±0.42 


0.9±0.4 


-19.95±0.42 


0.07±0.03 


A9.1 


-19.97±0.41 


0.5±0.2 


-23.81±0.41 


2.4±1.0 


A31.1 


-19.76±0.19 


0.4±0.1 


-17.62±0.19 


0.008±0.001 


M0304 


-22.46±0.20 


5.0±1.0 


-23.44±0.20 


1.7±0.3 


M0359 


-19.82±0.35 


0.8±0.3 


-21.42±0.35 


0.3±0.1 


M0416 


-19.65±0.19 


0.4±0.1 


-20.60±0.19 


0.1±0.02 


M1115 


-13.65±0.40 


0.004±0.002 


-16.43±0.40 


0.004±0.002 


Ml 206 


-21.55±0.16 


1.7±0.3 


-21.50±0.16 


0.24±0.04 


M1311 


-18.22±0.72 


0.2±0.1 


-19.95±0.72 


0.08±0.06 


M2031 


-23.39±0.45 


9.3±4.2 


-23.00±0.45 


1.0±0.4 


M2129 


-13.28±0.29 


0.002±0.001 


-22.07±0.29 


0.8±0.2 



sation species, but we do not detect any offset in velocity between 
the lines. 

To determine if velocity offsets are seen in our sample, we list 
in Table|16|the velocity offsets derived from the difference between 
the average emission and absorption redshifts. When the velocity 
offset is positive, one can interpret this as an outflow. Within the 
errors, there is no evidence for large velocity offsets in our data. 
A more detailed study would require higher S/N data in particu- 
lar to study the widths of the ISM lines. As the outflows in LBGs 
are driven by feedback from star formation, the lower SFRs in our 
galaxies relative to LBGs could result in smaller velocity offsets. In 
contrast, in low-mass galaxies which have weaker potential fields 
even small SFRs can drive a stronger wind. In such cases, the spe- 
cific SFR is expe cted to correlate w ith the wind strengths as seen in 
2 ~ 1 galaxies (Kor nei et alj201 j) . 




0.00 0.05 0.10 0.15 0.20 
E(B-V) 



5.2 Gas phase reddening 

For five sources, more than one of the Balmer lines are detected and 
their flux ratios can be used to derive the gas phase reddening. As 
the ratios of emission lines are not affected by lens magnification 
and slit losses, we do not include those parameters in the calcula- 
tions. The expected emission line ratio in the absence of reddening 
is a function of the gas temperatures and densities. We adopt the 
values tabulated for te mperatures T = 10 00 0-20 000 K and den- 
sities n ~ 100 cm -3 in Brockle hurstl d 197 ll) . The electron density 
of the gas in the lensed galaxies can be determined from the flux 
ratio of the [O II] AA3727,3730 doublet. The ratios of the lines are 
consistent with a low electron density <300 cm -3 in all cases. 

To derive the gas phase reddening, E(B — V / ) gas , we use all 
of the Balmer line pairs available for each of the galaxies and use 
a weighted average to determine the reddening. In practise this im- 
plies that the reddening is determined by the Balmer line ratio with 
the smallest uncerta inty. We assume a starburst extinction curve 
dCalzetti et alJl200oT) . and list the resulting reddening in Table [T6l 
The reddening derived for M2031 is very uncertain, but consistent 
with zero which was also found in the SED analysis in Sect l4.ll 

In the following analyses, all the emission lines fluxes in Ta- 
bles 5-13 are corrected for the intrinsic reddening E(B — V) sas , 
while for those galaxies where we do not constrain this parameter 
from emission line ratios we use instead the reddening measured 
for the stars in Tablell4l and assume a conversion to t he reddening 
in the gas phase E(B- V) sas = E (B-V) stars /0.4A JCalzetti et all 
2000). Figure [6] shows the reddening in the gas phase as a func- 



Figure 6. Reddening in the gas phase as a function of the stellar reddening 
determined from the SED fi tting. The dashed lin e shows E(B — V) gs , B = 
E(B - y) stars /0.44 from lCalzettietalJfeOOOh . 

tion of the reddening of the stars determined from SE D fitting. The 
dashed line represents the scaling factor of 0.44 from lCalzetti et al] 
(2000), which appears to be valid for the lensed galaxies as well. 



5.3 Starburst or AGN ionisation 

The ionising radiation from either active galactic nuclei (AGN) or 
recent star formation gives rise to strong emission lines in galaxies. 
Whether the emission is caused by an AGN rather than strong star- 
bursts is usually investigated through emission line ratios, which 
are sensitive to the hard ionisation from AGN. Conventionally, 
when only the strongest rest-frame optical emission lines can be 
detected in hig h-redshift galaxies, the line ratios [O IIl]/H/3 ver- 
sus [Nll]/Ha dKewlev & Dopitall2002h can be used. However, in 
our sample all these strong lines are only detected in M0304 and 
M0359. Both these galaxies have line ratios characteristic of star 
forming galaxies. 

In M2031, the C IV AA1548,1550 lines are not detected to a 
level of 3 x 10~ 18 erg cm" 2 s"\ so the ratio C IV/Cm]<0.15 
indicates a softer ionising spectrum, whil e an AGN would res ult in 
a fraction of 2 as argued for the Lynx arc dBinette et alj2003l) . The 
line ratio for A31.1 is C IV/[ClII,ClIl] = 0.65±0.71 indicating a 
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higher ionisation parameter. In AGNs, the lines N V AA1238,1242 
and Hell A1640 are expected to be strong, and since these lines 
are not detected in any of your spectra it suggests that the main 
contributor to the ionisation is massive stars. 

With the numerous emission lines detected from the sources 
it is possible to derive the ionisation parameter, U, for the galax- 
ies. We use the [Om] AA4959,5007/[O II] AA3727,3730 combined 
with the oxygen a bundances to calculate l og U using the itera- 
tive calculations in lKewlev & Do pita (2002). The resulting values 
lie in the range between -3 and -2, which is higher than mea- 
sured in l ocal galaxies, but typic al for high redshift galaxies, both 
unlensed dirt) et alJ|2006l. l201Ch and lensed ones dHainline et al.l 
120091 iRichard et alj201 lahlRigbv et al.ll201 lb . 

5.4 Oxygen abundance 

Oxygen abundances can be determined from the strong rest-frame 
optical emission lines. These values are not magnification depen- 
dent, because the parameters are derived based on line flux ratios 
rather than absolute values. A more accurate method to determine 
abundances rely on the detection of temperature-sensitive emission 
lines which are typically faint and very rarely detected in high- 
redshift galaxies. In three of the lensed galaxies at 2 < z < 3.5, 
the temperature-sensitive ratios of [O III] A5007 and the Oin] 
AA1661,1666 doublet allows us to use th e direct T e methods, which 
is presented in lChristensen et all d2012l) . 

When the temperature-sensitive lines are absent from the 
spectra, one has to rely on the relations between strong emis- 
sion line ratios a nd the oxygen ab undances calibrated from di- 
rect T e methods dOsterbrockl 1 1 989h . The ratio defined as R23 = 
(fOII1AA3727,37 30 + [OIII]AA4959, 5007)/H/3 introduced by 
IPagel et al.l ( fl97^ is one of the most used calibrations, since the 
relevant emission lines are detectable from ground-based data out 
to z — 3.8. However, the R23 relation has two possible solutions 
with a high-metallicity and a low-metallicity branch. A wealth of 
authors have published separate calibrations of t he R23 diagnos- 
tics based either on photo-ionisation models (e. g. McGaugh 1991; 
iKobulnickv etafll 19991 : iKewlev & Dopitai r2002). or direct calibra- 
tions derived from temperature-sensitive lines (e.g. lAlloin et al.l 
ll979l : |Pilvugin & Thuarj 20051 : iNagao et all2006h . Typically, these 
calibrations posses an intrinsic scatter of ~0.2 dex, and in absolute 
terms vario us calibrations may diff er by up to 0.7 dex relative to 
each other dKewlev & Ellisori2008h . 

A degenerate solution can be avoided if other emission lines 
such a s [Nil] A6586 or [Neiii] A3869 are detected foenicolo et ail 
20021: IKewlev & Dopital I2OO2I : iPettini & Pagell 120041: iLiang et all 



2006). Also the [O m]5007/[O II]3727 ratio can be used to dis- 



criminate between the upper and lower branch of the calibration 
(Naga oetalJl2006h . We detect [Nil] A6586 in two of the spectra, 
while [Ne III] A3869 is detected in three galaxies. For seven sources 
we can determine the R23 ratio, including those wit h direct oxygen 
abundance determinations Ichristensen et al .112012b . 

For the A4. 1 source, the H7 flux is scaled to determine H/3. Its 
Ha to H7 ratio suggests a reddening of E(B — V) = 0.24±0.13, 
so correcting Ha for the intrinsic extinction, and given that H/3 is 
a factor of 2.86 fainter than Ha we derive a corrected H/3 flux of 
(37.6 ± 3.2) x 10 -18 erg cm -2 s _1 . The other emission lines are 
also corrected for this this reddening value. 

For some galaxies, [O III] A4959 is not detected, noisier be- 
cause of telluric absorption, or close to a sky emission line. In those 
cases, we use the relation J5 007 = 2.98 x 74959 t o determine 
the total flux of the doublet dStorev&Zeippenll2000l) . To derive 



the oxygen abundance, we use the ex pression for the upper an d 
lower branch of the R23 calibration in IPilvugin & Thuanl d2005h . 
The R23, 12+log(0/H) values, and the propagated errors are listed 
in Table QJo] For completeness both the upper and lower branch 
abundances are listed. Gas phase abundanc es are derived relativ e 
to the solar value of 12 + log(0/H) = 8.69 dAsplund et aT]|2009h . 
We can use the additional detected lines for some of the galaxies to 
constrain which of the uppe r or lower branch o f the abundance cal- 
ibration is valid. Following ICresci et all d2012h . we assume the up- 
per R23 branch is valid for galaxies with log([0 IIl]/[0 II])< 0.45 
and log([Ne IIl]/[0 II])< 0.6. Whether the upper or lower branch is 
preferred is indicated by the oxygen abundances listed in bold face 
numbers in Tablell6l 

To further justify the choice of the upper or lower branch of the 
metallicity calibrations, Table |16|also lists the oxy gen abundances 
using the Q3N2 calibration dPettini & P agel 2004) and Ne302 cal- 
ibration dNagao et alj|2006t) . In all cases, these calibrations agree 
with the preferred branch of the R23 relation. The various oxy- 
gen abundance determinations vary within 0.3 dex, which illus- 
trates xy^fHheJjnownoffs^t^ between various strong line calibra- 
tions dKewlev & Ellisonl2008h . 

In the case of M2031, M0304 and Abell 1689 ID 31 1, al l 
have direct abundances measurements (see lChristensen et alj2012f) 
consistent wi th those inferred from th e lower branch of the R23 
calibration in IPilvugin & Thuanl d2005h . 



5.5 Star formation rates 

We determine the SFRs using three methods: [Oil] and Ha emis- 
sion lines, and the UV continuum flux at rest-frame 1500-2800 A. 
The measured emission line fluxes are corrected for the slit losses 
and magnification factors (Table |3j, and we also correct the fluxes 
for the intrinsic reddening. From the sum of the flux in the doublet 
[Oil] A3727,3730 we calculate the extinction-corrected luminos- 
ity (I/[on],cor based on the gas phase reddening, E(B — ^) gaa ), 
and use the relation in lKennicuttl ( fT998T) to derive the intrinsic SFR. 
We divide the conversion factor by 1.7 to correct the SFR from a 
Salpeter IMF to a Chabrier IMF. In a similar manner we calculate 
the SFR based on the Ha lines, which is detected for three of the 
galaxies. Again the SFR conversion is divided by a factor of 1.7. 
The resulting values are listed in Tablell6l In addition, we also cal- 
culate the SFR from M2031 assuming that Ha is 2.86 times the H/3 
flux in case of a zero reddening. In the following section|6] we use 
the SFRs based on the [O II] AA3727,3730 doublet. Once corrected 
for reddening, S FRs based on [O II] and Ha are generally in good 
agreement (e.g. IKewlev et ai]|2004h . and Ha is known as a good 
tracer of the current ongoing SFR. 

As an alternative, we also calculate the SFR from the rest- 
frame UV lumi nosity at 1500-2800 A using the calibration in 
iKennicuttl {l~998), again correcting the scaling relation for the IMF 
difference by dividing by a factor 1.7. We make use of the best fit 
spectral models available for all of the galaxies in our sample, cor- 
rect these for the reddening E(B — V%tars, and determine the UV 
flux at the rest-frame 1500 A. The resulting SFRs are listed in Table 
M 

The uncertainties for the SFRs are propagated and contain the 
uncertainties measured either for the [O II] flux or the UV contin- 
uum as well as the uncertainty for the slit loss and magnification 
factors. When two or more estimates of the SFRs can be derived, 
we find that the values are consistent to within a factor of two. In 
contrast to the SFRs, the specific SFRs, i.e. the SFRs divided by the 
stellar mass, do not depend on the magnification factor, slit losses 
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Table 16. Physical properties of sources with emission lines 



Arc 


A4.1 


A31.1 


M0304 


M0359 


M0416 


M1311 


M2031 


Section 


• 

Velocity dispersion a (km s ) 


17±1 


22±2 


56±1 


53±2 


65 ±9 


53±3 


43±1 


kill 
|5.1.1| 


A V (2cm - 2 a bs) (km S _1 ) 




95±131 


-10±55 








80±90 


|5.1.2| 


SFR([OII]) (M yr- 1 ) 


0.7±0.1 


1.8±0.3 


55±11 


8.1±2.7 


4.7±0.7 


3.0±1.8 


7.1±3.3 


EH 


SFR(Ffa) 


0.4±0.1 




16±3 


1.5±0.5 






14.5±6.5 




E(B - V0 gas (mag) 


0.24±0.13 


0.37±0.28 


0.23±0.01 


0.29±0.01 






0.03±0.33 


H3 


log U 


-2.7 


-2.1 


-2.8 


-3.1 


-2.4 


-3.1 


-2.1 


m 


log i?23 


0.81±0.04 


0.85±0.05 


0.94±0.001 


0.67±0.01 


0.81±0.04 


0.19±0.06 


0.91±0.02 




12+log(0/H) (upper branch) 


8.33±0.08 


8.45±0.04 


8.22±0.01 


8.18±0.05 


8.44±0.06 


8.66±0.14 


8.39±0.01 


EH 


12+log(0/H) (lower branch) 


7.91±0.11 


7.63±0.10 


8.16±0.01 


7.48±0.04 


7.73±0.11 


7.13±0.01 


7.74±0.03 




Z/Zq (lower; upper) 


0.2; 0.4 


0.1; 0.6 


0.3; 0.6 


0.3; 0.3 


0.1; 0.3 


0.03; 0.9 


0.1; 0.5 




12+log(0/H) (03N2) 






8.04±0.01 


8.48±0.01 










12+log(0/H) (Ne302) 


8.24±0.07 


7.85±0.15 


8.20±0.01 








7.56±0.11 





Notes. 

The listed SFRs assume a Chabrier IMF as described in the text. For completeness we list both the upper and the lower branch calibrations for the R23 
diagnostics. The values in bold face represent the preferred values based on the ratios of either [Ne III]/[0 II] or [O III]/[0 II]. 



or choice of IMF. The inverse of the SSFRs can be taken as the 
formation time-scale of the galaxy given that the SFR is constant. 
For a galaxy with a small SFR and high stellar masses, the corre- 
spondingly low SSFR suggests that the SFR must have been higher 
in the past. In particular, M2129, the oldest galaxy in our sample, 
shows evidence of only little present ongoing star formation. The 
[O II] AA3727,3730 doublet is detected at a level corresponding to 
a SFR = 0.8 Mq yr -1 , implying that much stronger star forma- 
tion has taken place in the past to build up the stellar mass in this 
massive galaxy. 



6 FUNDAMENTAL RELATION FOR STAR-FORMING 
GALAXIES 

6.1 Mass-metallicity relations 

Observations of mass-metallicity relations have revealed an evo- 
luti on with increasin g redshifts from local gala xies in SPS S 



1 Tremonti et al. 2004), fiel d galaxies at z = 0. 7 ( Savaalio et al 
2005), and out to z > 3 jMaiolino e t al. 20081: iMannucci et al 



2009) . iRichard et al] feoi 11 an " analyse a larger sample of lensed 
galaxies at 1.5 < z < 3.3, and find the galaxy metallicities at 
z ~ 3 are slightly larger than field galaxies. 

Figure |7J shows that the lensed galaxies exhibit a huge scat- 
ter compared to the established relations, and we do not find any 
clear correlation. The highest stellar mass galaxy, M0304, belongs 
to a merging system, and therefore naturally has a lower observed 
metallicity relative to the expected value from extrapolation of the 
mass-metallicity relations. 

Since the mass-metallicity relations are based on various emis- 
sion li ne diagnostics, which e xhibit offsets with respect to each 
other ( Kewley & Ellison 2008), it is not straight forward to analyse 
the evolution of the mass-metallicity relation with redshift when 
different galaxy samples are involved. Such offsets in abundance 
diagnostics can partly explain the reason why the lensed galaxies 
do not follow any of the correlations at a given redshift. 




8 9 10 

log (M./M s ) 



Figure 7. Oxygen abundances of the lensed galaxies as a function of the 
stellar masses, corrected for lens magnifications. The black circles show 
the oxygen abundances from the preferred branch of the R23, relation in- 
cluding the direct measurements for A3 1.1, M0304 and M203 1 , while the 
grey error-bars show the location of the galaxies if the alternate branch 
of the .R23 relation is adopted for all galaxies. The mass-metallicity rela- 
tion for Lyman break galaxies at 2 = 3 — 4 is shown by the solid line 
iMaiolino eta l. 2008; M annucci et alj|2009l) . and galaxies at 2 = 0.7 by 
the dashed line jSavaglio et alj2005l) . 



6.2 Fundamental relation at low stellar masses 

The fundamental relation between galaxies SFRs, oxygen abun- 
dances, a nd stellar masses calibrated for low-redshift, h igh-mass 
galaxies dLara-Lopez et ai] |2010l : IMannucci et ail [2010) have re- 
vealed a larger scatter when extending relation to lower stel- 
lar masses. In particular, GRB host galaxies with M* <; 
10 9 Mq s how a larger scatter relative to the funda men- 
tal relation jKocevski & Westl l201ll ; IMannucci eTal] 1201 ll) . and 
IMannucci et al. ( 201 1 ) provide an extension of the fundamental re- 
lation at lower stellar masses. Recent observations show that also a 
low stellar m ass supernova Type la h ost galaxy at z = 1.55 follows 
this relation ( Frederiksen et al. 2 0121) . 
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Figure 8. Offsets of the derived oxygen abu ndance compared to th e pre- 
dicted one from the fundamental relation in Mannucci et al] [20\ ll) . The 
red circles show the location of the preferred R23 branch with the calibra- 
tion in lMaiolino et al. 1 2008), a nd using the dire ct measurements for A3 1.1, 
M0304 and M2031 in lChristensen et alj 120121) . A location above zero im- 
plies that the measured metallicity is higher than expected from the relation. 
The grey squares show the offsets of 13 lensed g alaxies at z = 2 — 3 for 
which either the R23 or the N2 calibration from Maiolino et al. (20(3) were 
used jRichardetallteOllall The grey triangle s show offsets of 10 lensed 
galaxies at 2 = 1 — 2 in lwuvts et all J2012bl) for which the oxygen abun- 
dance was derived from the N2 calibration in Maiolino et al. (2008). The 
outlined diamonds represent 1 8 GRB hosts at z < 1 {Mannucci et al j201lh . 
and the star symbols represent unlensed and more massive z ~ 3.5 galaxies 
I Maio lino et alj2008h . 



in the literature follow the fundamental relation when their metal- 
licitie s are 12+log(Q/H) > 8.0, and stellar masses larger tha n 10 10 
M dDessauges-Zavadskv et~ail201 ltlRichard et al]|201 id) , while 
lensed galaxies with stellar masses below ~ 10 10 Mq have in ma- 
jorit y larger metallicity than expected from th e fundamental rela- 
tion dRichard et al ] |2011al . lWuvtsetalJl2012bh consistent with our 
findings. 

If we assume that the oxygen abundance in the fundamental 
relation equation depends only on the SFR, the galaxies would fall 
within 0.5 dex of a relation. Similarly, a smaller scatter with respect 
to the fundamental relation can be obtained if the SFR is not in- 
cluded, and only the stellar mass is a dependent variable. However, 
the minimum scatter with respect to the relation is maintained with 
a combination of depe ndence of both the stellar mass and the SFR. 
Mannu cci et al. lime]) defines a projection of the fundamental re- 
lation: n a — log(M») — alog(SFR), where a = 0.32 gives the 
smallest scatter. For the lensed galaxies with stellar masses below 
10 10 Mq, we determine the smallest scatter at a = 0.36. 

Since the sample of lensed galaxies is so small compared to 
the 140 000 SDSS galaxies, we can only suggest an alternative fit 
for the fundamental relation valid for low mass galaxies. The oxy- 
gen abundances for all the le nsed galaxies in the rig ht panel of 
Fig. [Hand the GRB hosts from lMannucci et aU d201 1 ) are fit with 
the fu nctional form of the fundamental relation in 
(2010). The resulting fit is: 



Mannucci et al.l 



12 + log(0/H) = 8.52 + 0.26m + 0.013s + 0.025m 2 
-0.026ms + 0.013s 2 , 



(1) 



If we calculate the differences bet ween the oxygen abun- 
dances in Table [16] and the calibration in Man nucci et al.l d201lh . 
we find that the galaxies exhibit a large scatter compared to the re- 
lation. A potential risk when analysing and comparing abundances, 
is that there is a significant offset b etween various emission line 
diagnostics dKewlev & Ellisonl|2008l) . We t herefore also derive th e 
oxygen abundances using the calibrations in lMaiolino et a 
which correspond to th e abundances used in t he derivation of the 
fundamental relation in lMannucci et al.l d20ld) . The offsets are il- 
lustrated in Figure [8] Other low -mass lense d galaxies at z = 1 — 3 
dRichard et alTl201ld : | Wuvts et al. 2012b) illustrated by squares 
and triangles in Fig. [8] show a similar trend of not following the 
fundamental relation. The standard deviation of the offsets from the 
relation is 1 dex, and the galaxy with the lowest stellar mass, A4.1, 
is offset by ~1 dex. In any case, the high abundance of this galaxy 
is unexpected given its low stellar masses. The basis of the selec- 
tion of the upper branch for the R23 abundance calibration is that 
the flux rat ios of [Ne III]/[Q I I ] and [O III]/[Q II] suggest a higher 
abundance dNagao et al . 2006; Cr esci et alj|2012l) . 

In addition to the large scatter, the locations of the galaxies 
have a slope of -0.3 per dex in log M*, which is statistically sig- 
nificant with a Spearman nonparametric correlation test giving a 
probability of 0.98. Th e calibration of the fundamental relation in 
lLara-Lopez et al] d2010l) also produces a slope and a large scatter. 

Our finding of a discrepant location of the lensed galaxies 
compared to the funda mental relation for star-forming galaxies in 



Mannucci et al. (2010) and the extension to low-mass galaxies in 
Mannucci et al] ((2011) does not imply that a fundamental relation 
is not present. It merely suggests that the scatter is significant when 
low-mass galaxies are analysed. Massive lensed galaxies analysed 



where s = log(SFR) and m = log(M„) - 10 in solar units. With this 
alternative calibration, the scatter is 0.3 dex and the slope is absent 
as demonstrated in Fig. [9] A Spearman test gives a probability of 
0.3 for a correlation. 

One concern is that the R23 calibration is not an accurate 
metallicity t racer given the large offse ts between various strong line 
calibrations dKewlev & Ellisonl2008l) . A comparison of the oxygen 
abundances from R23 and the direct T e method for metal-rich, lo- 
cal SDSS galaxies suggests that the T e -based log(0/H) abundances 
are lowe r by up to 0.6 dex c ompared to the upper branch of the R23 
relation dLiang et al.lr2007t) . However, f or the galaxies in our sam- 
ple we find that the R23 calibration in Pilvug in~& Thuanl (2005) 
is in agr eement with direct T* e b ased abundances for three of the 
galaxies dChristensen et all l2012h . For the other four galaxies we 
can only assume that this is valid too. 

Since the lens magnification as well as the slit loss correction 
affect both the SFR and stellar masses equally, any uncertainty in 
the magnification factor and correction factor has little impact on 
the location of the galaxies relative to the fundamental plane. The 
only assumption is that the emission lines are uniformly distributed 
over the entire galaxy, such that the emission line flux follows 
the continuum flux. Metallicity gradients are known to be present 
in disk galaxies, so we likely overestimate the integrated oxygen 
abundance slightly. Considering that gradients are typically a few 
times 0.1 dex over the face of a large (local) galaxy, and that the 
lensed galaxies in this study are compact sources dominated by at 
most a few luminous H II regions, any metallicity gradient cannot 
explain an integrated 1 dex offset relative to the fundamental rela- 
tion. 
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Figure 9. Offset in measured abundance compared to the new calibration 
of the fundamental relation given by Eq. (1). The symbols present the same 
data as in Fig. [8] From the entire sample of lensed galaxies combined with 
low redshift GRB hosts, we measure a standard deviation of <r = 0.3 dex. 



7 DISCUSSION AND SUMMARY 

We have presented spectroscopic data and determined redshifts for 
12 strongly lensed galaxies behind various massive galaxy clus- 
ters. Three of the redshifts were previously determined by other 
studies. Although the selection of the lensed galaxies is based on 
optical data, and therefore makes a preferential selection of UV 
bright high-redshift galaxies, we still achieve a sample of various 
galaxy properties (young and old, reddened and non-reddened, star- 
forming and non-star-forming), because we do not include colour 
selections. The lensing magnification effect allows us to carry out 
detailed investigations of galaxies that are on the average 1-2 or- 
ders of magnitude fainter and less massive than galaxies studied 
in unlensed field populations. The galaxies are also intrinsically 
fainter than other lensed galaxies studied in the literature, sim- 
ply because we targeted galaxies with observed magnitudes ~24, 
which is fainter than the ones investigated previously. 

Using the stellar masses, SFRs and oxygen abundances, we 
find a breakdown of the fundamental relation for star-forming 
galaxies at very low stellar masses (M < 1U 1O M ) with offsets 
of up to 1-2 dex, w hich implies that one cannot use the calibration 
in iMannucci et al] bOld) to derive any missing parameter (SFR, 
M», or oxygen abundance) when two other parameters are known 
for high-redshift galaxies. Contrary to the suggestion that the rela- 
tion is not valid at z > 2.5, we find one galaxy at z = 3.5 with a 
metallicity only 0.4±0.1 dex too low compared to the expectation 
from the fundamental relation. We provide an alternative calibra- 
tion extending the low-mass end of the fundamental relation. 

Our findings of the SFRs from various indicators (UV con- 
tinuum versus either [O II] or Ha emission lines) agree with each 
other within a factor of two after correcting for the intrinsic redden- 
ing. Such spreads in the SFR estim ates from v arious diagnostics are 
known for LBGs fcrb et alj[2003h . However. Iwiivts et al] J20123) 
find that the SFRs after dust correction in four lensed z = 2 galax- 
ies agree under the assumption that the extinction la w is steeper 
than a Calzetti law as observed for unlensed LBGs (Red dv et al] 
2010). In the lensed galaxies studied here, the UV-based SFRs 
are on the average lower than inferred from emission lines, which 
would support the suggestions of a steeper extinction curve, such 



that the UV flux should be corrected by a larger factor to match the 
SFRs based on emission lines. However, even if the SFRs are over- 
estimated by a factor of two, this would make no impact on the off- 
set of the lensed galaxies with respect to the fundamental relations 
between the SFR, stellar mass and oxygen abundances explored in 
Sect.[6l 

We find no evidence for velocity offsets between rest-frame 
UV and optical emission lines and UV absorption lines when both 
are detected, since the offsets we derive are only significant on the 
la level. Observations of high-redshift galaxies have interpreted 
offsets between Lya and UV absorption lines or optical emission 
lines as indication of galaxy outflows of up to several 100 km s , 
even though the velocity offset of Lya emission lines can be ex- 
plained by radiation transfer effect in expanding neutral gas without 
the need for extreme velocities. The combination of accurate red- 
shift determination from intermediate resolution spectra combined 
with a full spectral coverage is essential to investigate individual 
galaxies for these kinematic effects. High-velocity outflows from 
galaxies can also be inferred from the extension in velocity space 
of the ISM absorption lines, but in the low S/N ratio spectra investi- 
gated here, we cannot determine if the absorption lines span a large 
velocity range. 

One of the driving questions behind this investigation is which 
galaxies have strong emission lines that are detected with shallow 
spectra. Galaxies without bright knots (or equivalent to H II re- 
gions), do not have strong emission lines. A8. 1 and Ml 1 15 show no 
emission lines at all, while M1311 has only very faint lines. Abell 
1689 ID 8.1 is a giant arc, which has a blue continuum suggest- 
ing a young galaxy with ongoing star formation. With a redshift 
z — 2.68, it should have detectable emission lines, also beyond 
Lya. The Lya line itself could easily be absorbed as it is only de- 
tected in about 25% of LBGs. On the other hand, other lines could 
be present but just absorbed or obscured. As seen in the M1206 arc, 
we only detect a single but very strong emission line from [O III] 
A5007, while all the other strong emission lines in the large wave- 
length range spanned by X-shooter fall behind strong sky emission 
lines, or absorbed by telluric lines. In the case of Ml 1 15, which is 
a low-surface brightness arc, the photometric redshift is either low 
or high, so emission lines are expected. However, the SED of the 
galaxy suggests an old stellar population, where emission lines in 
H II regions are no longer produced. 

Alternatively, we can compare the relative luminosities at the 
rest-frame UV and optical. Not surprisingly, the galaxies with low 
absolute UV luminosity relative to that in the rest-frame (Ml 115 
and M2129) do not exhibit strong emission lines. When comparing 
the rest-frame UV and optical luminosities in Table|15l most galax- 
ies can be separated into UV bright or UV faint relative to their op- 
tical luminosities. We find that the galaxies without emission lines 
are those with fractional luminosities L/L*(UV)< L/L* (optical). 
The only exception is again A8.1. 

In future studies when large numbers of gravitational lensed 
candidates are detected in multi-band images covering large areas 
on the sky, we intend to focus not only on galaxies with obvious 
clumpy morphologies indicating the presence of H II regions and 
on-going star formation, but also on the galaxies with steep spectral 
slopes. 
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